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PREFACE 


It is now well recognized that certain types of macromolecules, 
when added in relatively small concentrations to solutions flowing in 
turbulent motion, possess the property of reducing the friction drag 
of the fluid in which they are suspended or dissolved. As a result, 
the addition of these substances can significantly increase the capacity 
of a system to transmit fluids at normally turbulent flow rates. 

The feasibility of utilizing this effect to reduce the time 
required to perform certain critical refueling operations encountered 
in military logistics support prompted this investigation. That the 
tactical vulnerability of a combatant force is at its peak during 
periods of logistic replenishment is well recognized by military logis- 
ticians. This is particularly true during periods involving underway 
refueling operations from a mobile replenishment group at sea, midair 
refueling of aircraft and the refueling of aircraft aboard aircraft- 
carriers at sea. 

The original scope of this investigation envisioned an examina- 
tion of the drag reducing effects of a number of viscoelastic polymers 
in several of the more common military fuels in use. It was further 
recognized that the effect these additives would have on the quality of 
mie tuel would be of paramount interest. However, because of Ehe time 
required to design, assemble and calibrate a flow system which would be 
Suitable for investigation of the drag reduction phenomena, this scope 
was limited to the examination of a single high molecular weight poly- 
isobutylene (Vistanex L-200) added to military jet fuel (JP-4). The 
effect of this additive on fuel quality was not investigated. 

A variable pressure, single pass flow system was designed and 
constructed for use in this study. This type of system was selected in 
order to maintain similarity with tactical refueling operations. 
Further, the effects of mechanical degradation could best be observed 
in a Single pass system since the fluid makes only one pass through the 
test section. Finally, flow behavior over a wide range of the turbulent 


flow regime could be observed yet fluid consumption and equipment cost 


cat 





ia 
could be minimized. The apparatus was found to perform well in both 
the laminar and turbulent flow regimes when tested with liquids of well 
known fluid properties. 

All drag reduction experiments were accomplished by observing 
the flow behavior of solutions containing 50 wppm of the high molecular 
weight Vistanex L-200 in military jet fuel. The results of the study 
indicate that a marked reduction in the viscous drag of jet fuel Flow- 
ing in the turbulent regime can be achieved. It is anticipated that 
the time required to transmit fuel could be better than halved. 

The polymer tested was found to be extremely sensitive to 
mechanical degradation. While the results of this study are not conclu- 
Sive, a striking similarity between what has been termed "the diameter 
effect" and the effects of mechanical degradation was observed. A cor- 
relation suitable for dynamic scale-up to full sized flow systems was 
frustrated by this degradation/diameter effect. Further, it was found 
that mechanical devices producing extremely high localized shear rates 
could completely destroy the drag reducing effectiveness of the additive. 
Therefore, it was concluded that this type of a device should be avoided 
in any full scale flow system or, alternatively, injection of the 
additive should be accomplished downstream from such a device. 

There were strong indications that a suppression of turbulent 
fluctuations might well be the main characteristic of drag reduction. 
Additional work with more precise equipment was recommended which would 


serve to verify this tentative conclusion. 
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CHAPTER 


INTRODUCTION 


The objective of this study was to investigate the effects on 
fluid flow caused by the addition of a minute quantity of a high 
molecular weight polymer to military jet fuel (JP4). The flow was 
observed in the turbulent flow regime corresponding to a range of 
Reynolds numbers from lo” to 10°. In particular this thesis is con- 
cerned with the characteristics of frictional drag (or its reduction) as 
influenced by such an additive under conditions of flow in cylindrical 
tubes. The primary emphasis of this study has been to seek the factors 
which tend to influence the effectiveness of such additives in a common 
fuel in actual use, rather than to search for the mechanisms which might 
be responsible for the cause of drag reduction. Further, it was recog- 
nized that time would permit only a very limited examination of the 
influence of these additives since it would first be necessary to design, 
assemble and calibrate a flow system suitable for the investigation of 
the drag reduction phenomenon. However, it was hoped that this effort 
might become the basis for future work in this area since a calibrated 


system would be readily available. 


LITERATURE SURVEY AND DISCUSSION 

Cayerar Reducer onion ne dm luc terme dircde er edMcrlen ser 
"fluid friction reduction" has been used to describe a rather broad 
Spectrum of rheological phenomena since first introduced in the technical 
literature in 1949. It is therefore necessary to consider two quite 
distinct classifications of fluids, purely viscous and viscoelastic, with 
which the term is frequently associated in order to clarify the defini- 
tion as used in this thesis. The first classification consisting of 
fluids commonly described as "purely viscous" can be further subdivided 
Mato two categories. The fluids of the first category are said to obey 
Newton's Law of Viscosity and are broadly classified as "Newtonian 
Fluids." The behavior of such fluids can be described by the general 


power law or "Ostwald-de Waele" model (3): 





2 
tT=-m (dv /dy)” Pil 


where n = 1 and m is a constant andm = u 

This relationship states that the shear force per unit area (shear 
stress) varies directly with the negative of the velocity gradient. A 
more useful relationship, since it describes the shear stress/shear 
rate characteristics of laminar flow in cylindrical tubes in terms of 


measureable or derived quantities, is: 
Cs DAP/4L = yp (8V/D) [1-2] 


Or 


RAP/2L = wu (4Q/nR°) [1-3] 


Results of extensive experimental tests of these fluids in turbulent 
flow through sand-roughed pipes were presented by Nikuradse in 1933. 
Friction factor data are conventionally presented as functions of 
Reynolds numbers in dimensionless plots using logarithmic scales with 
relative roughness as a secondary variable. Perhaps the most well 
known plots of this type are the "Moody Diagrams" (14) found in common 
engineering handbooks. (9) 

Fluids in the second category, "purely viscous," are classified 
as "non-Newtonian" and can generally be described by the power law rela- 
tion, equation [1-1], but without the restrictions that n = 1 and m be 
constant. Similarly, in terms of physical quantities more amenable to 


direct measurement: 


eon ieee [2(3 + 1/n) V/D]" oe 


After measuring AP as a function of V (or flow rate from which an aver- 
age velocity can be derived) a plot of the variables DAP/4L versus 

8V/D on log-log paper can be constructed which will reveal to which of 
these two viscous categories a fluid belongs. The popularity of the 
Ostwald-de Waele power law as a means of describing viscous fluid 
behavior now becomes readily apparent since it has been shown to cor- 
rectly portray the relationship between shear stress and shear rate for 


many fluids over the broad range of shear rates commonly encountered. 





Extensions of the experiments conducted by Nikuradse which 
include the behavior of non-Newtonian fluids have been the subject of 
more recent investigators. (4, 22) By using a "generalized" Reynolds 
number which has been shown to be a function of power law coefficients, 
m and n, Dodge and Metzner (5) have revealed the nature of the friction 
factor, "generalized" Reynolds number relationship for flow of non- 
Newtonian fluids through smooth cylindrical pipes. The results of the 
Dodge and Metzner studies indicate that: 

(a) The generalized Reynolds number at which transition 
from laminar flow to turbulent flow is observed increases slowly from a 
value of 2100 for Newtonian fluids as the flow behavior index, n, 
decreases. (n is customarily defined as the index of non-Newtonian 
behavior. ) 

(b) At any given generalized Reynolds number in the turbu- 
lent flow regime there will be a decrease in the friction factor 


observed corresponding to a decrease in the flow behavior index, n. 


These results might suggest behavior not totally unlike that associated 
with the drag reduction phenomenon. _ 

Since the friction factors for purely viscous non-Newtonian 
fluids usually will be lower than those of Newtonian materials if n<l, it 
would appear as if the additions of solids or polymeric materials to a 
Newtonian solvent to yield a non-Newtonian fluid might serve to reduce 
pressure drop under turbulent conditions. In fact, however, the thicken- 
ing action promoted by these additives generally will outweigh the fric- 
tion factor reduction so that the actual pressure drop required to 
develop a certain flow rate will almost always be increased. As correctly 
commented upon by Savins (21) and Fabula and Hoyt (6), it is important 
to recognize the difference between non-Newtonian behavior and the drag 
reduction phenomenon, which is to be discussed shortly. 

The high molecular weight, random-coiling polymers in solution 
appear to exhibit elastic characteristics in addition .to the viscous 
characteristics which describe the behavior of most fluids. Hence, the 
second broad classification of fluids to be considered is termed 
"viscoelastic." It might be well to recall at this point that the prop- 


erties of elasticity and viscosity are defined in terms of the response 





of a fluid to an applied stress. The elastic response describes the 
tendency of a fluid to recover its former configuration upon removal 

of an applied stress. This recovery would be complete were it not for 
the viscous flow which begins when the stress is applied. The viscous 
flow is never reversed and may continue once the stress has been removed 
owing to the inertial properties associated with the fluid. Thus, an 
essential distinction between elasticity and viscosity is that the 
former describes a reversible process while the latter describes an 
irreversible property. 

As might well be expected, very dilute solutions containing these 
sescoclastic polymers exhibit chandeterisrics very closely identisived 
with the fluid properties described as "Newtonian" by conventional vis- 
cometry i.e., n = 1, m= in Eq. [1-1]. Owing to the very low concen- 
trations of polymer additives, the measured density and viscosity of a 
solution are very nearly the same as those of the pure solvent. At: 
higher levels of concentration, they appear to impart non-Newtonian 
characteristics and while the assignment of "viscosity" as a fluid 
property to such fluids 1s somewhat nebulous, a thickening of the fluid 
is readily apparent. While investigating the characteristics of non- 
Newtonian fluids, both Dodge and Metzner (4, 5) and Shaver and 
Merrill (22) selected some additives which are now known to produce 
viscoelastic solutions. An apparent anomalous behavior was observed for 
the viscoelastic fluids when compared to the purely viscous fluids with 
approximately the same index of non-Newtonian behavior. 

It is at this point that authorities appear to disagree on the 
meaning of "drag reduction." Shin (23) contends that if a polymer solu- 
tion is "concentrated," it will display considerable gradient dependency 
of viscosity (the characteristic most commonly associated with non- 
Newtonian behavior) and that the assignment of viscosity for purposes 
of evaluating a Reynolds number is always equivocal. This suggests that 
the "flow anomaly" associated with the concentrated viscoelastic solu- 
tions observed by Shaver, Dodge, et al., might well stem from an imper- 
fect model used to evaluate viscosity and generalized Reynolds number 
rather than an actual reduction in friction stress. In contrast, the 


Newtonian characteristics of the solvent are retained when dilute solu- 





tions containing viscoelastic polymers are formed and the classical 
properties associated with viscosity and Reynolds number can be unequiv- 
ocally evaluated in a conventional sense. Further, in very dilute 
polymeric solutions, Jae fluid properties generally used to describe 
fluid flow are very nearly equal to the properties of the pure solvent 
when subjected to conventional measurement techniques. Thus, the depar- 
ture from conventional Flow behavior associated with these extremely 
dilute solutions can very clearly be labeled "anomalous." This is not 
to suggest that these same characteristics do not occur in concentrated 
viscoelastic solutions. Indeed, the observations of Dodge (4) and 
Shaver and Merrill (22) and more recently those of Savins (21) and 
Rodriguez, Zakin and Patterson (20) strongly suggest that frictional 
stress reducing characteristics separate and apart from non-Newtonian 
behavior of purely viscous fluids are present in these concentrated 
solutions. However, from a practical standpoint, the final criterion 
of the effectiveness of such an additive must be measured by its fric- 
tion reducing abilities as compared to the drag of the pure solvent in 
which it is dissolved. Thus, whether a final solution is considered 
Newtonian or non-Newtonian, dilute or concentrated, it is the frictional 
drag of that solution as compared to the frictional drag of the pure 
solvent which shall be used as the criterion to define drag reduction 


imecils thesis. 


(b) Causes and Effects of Drag Reduction Phenomena While 
there have been a number of attempts to isolate the precise mechanisms 
which give rise to the drag reduction phenomena noted in the literature, 
the exact mechanism is still not known. This, of course, is understand- 
able considering that a microscopic phenomenon is being interpreted on 
the basis of macroscopic data fitted to a macroscopic or microscopic 
model. 

The questions to be examined if one is to gain insight into the 
mechanisms which give rise to drag reduction, can now be stated more 
explicitly as: how the macroscopic fluid properties and flow behavior 


are related to the microscopic configuration and motion of the molecules. 


and, how the measurable macroscopic quantities such as fluid density, 





6 
viscosity and flow velocity can be shown to be consistent with a macro- 
scopic or microscopic model of fluid flow. 

The earliest published technical report available which attributes 
drag reduction to dilute polymer solutions appears to be that of B.A. 
Toms appearing in the Proceedings of the First International Rheological 
Congress. (25) Indeed, because of the pioneering research he conducted 
jmethis tield, 1 as often vefenpcdemorascsthes moms Ettect.|) a7 ee lnemes 
paper, Toms presented data which indicated an observed increase in the 
rate of flow of monochlorbenzene at a constant pressure gradient caused 
by the addition of polymethylmethacrylate. The remarkable feature of 
these results was that this increased flow rate was a phenomenon 
observed only in the nonlaminar flow region. 

Oldroyd (15) concurrently presented a paper suggesting that these 
unexpectedly high flow rates observed by Toms be explained in terms of 
anomalous behavior of a thin laminar layer at the tube wall when the 
mainstream of flow becomes turbulent. Oldroyd reasoned that, if the 
fluid is a solution containing a long chain linear polymer, an external 
constraint may be imposed by the wall of the flow system on the ways in 
which these long chain molecules can rotate. It might therefore be 
possible that an abnormally mobile pseudo-laminar sublayer, or boundary, 
of a thickness comparable with molecular dimensions could exist at the 
wall. This effect would then be macroscopically evident as an effective 
slippage at the wall. The flow rates as determined for conventional 
Newtonian fluids would have to be increased to account for this slippage. 

Although the proposed velocity of slip might well explain the 
macroscopte effect of drag reduction, a departure from the accepted 
model of fluid flow is implied if this theory is used to describe the 
microscopic cause. In contrast to the model used to describe ordinary 
flow of a Newtonian fluid, it was visualized that the long chain poly- 
mers moved with translational motion, hence would slide or slip along 
the side of the wall. Subsequent investigation however, has failed to 
reveal any translational slippage. To the contrary, flow experiments 
conducted in both the laminar and turbulent flow regimes appear to 
confirm a flow model based on a zero velocity at the wall. Despite the 


fact that Oldroyd's slip mechanism as a cause for drag reduction has not 





been given serious consideration for some time, there are relatively 
current texts which continue to refer to slip as being "important in 
some non-Newtonian flow problems." (3) 

A second concept which might explain the mechanism or cause of 
drag reduction was first advanced by Davies, Ward, Atkinson and Blair 
(17) in reaction to Oldroyd's boundary layer theory. These gentlemen 
suggested the possibility that the Toms Effect is tied to the elasticity 
of polymer solutions, and although variations of this theory appear to 
be the most widely accepted today, the exact mechanism by which the 
elastic properties influence drag reduction is still not fully ee 
stood. In general, this theory holds that the elastic properties of 
the fluid absorb the kinetic energy of turbulent motion, thereby reduc- 
ing the fluid transmission energy lost to turbulence. In other words, 
the irreversible friction losses due to turbulence are reduced because 
the energy associated with turbulent fluctuations has been reversibly 
stored in the fluid. 

A third mechanism, the formation of a viscosity gradient, was 
apparently originally proposed by Shaver and Merrill (22) though it is 
discussed in some detail by Shin (25) under the heading "Anisotropic 
Viscosity." This theory stems from Merrill's suggestion that the local 
viscosity normal to the direction of shear might well be substantially 
increased through elongation of the randomly coiled macromolecules even 
though the viscosity parallel to the direction of shear is very nearly 
identical to that of the solvent. Thus, the high perpendicular vis- 
cosity can reduce turbulent fluctuations in that direction. Since the 
turbulent fluctuations are randomly directed, the orientation of the 
elongated molecules will also be random. 

During the early 1960's, the Navy became particularly interested 
in the drag reduction phenomenon because of its possible influence on 
ship and torpedo performance. Fabula and Hoyt (6) of the Navy's 
Ordnance Test Station at China Lake, California, were particularly 
instrumental in surveying many water soluble drag reducing additives. 
As a result of their experiments conducted with aqueous polymer solu- 
tions, general rules were developed as to the type of material likely 


to be effective. In general, it was concluded that at least three 





Significant parameters affect the ability of a polymer to lower the 
turbulent frictional resistance of the fluid in which it is dissolved: 
linearity, molecular weight, and solubility. The "long chain" mate- 
rials have an essentially unbranched molecular structure of high 
molecular weight when added to solvents in which they can be fully 
dissolved, appear to produce the most marked drag reduction effect. 

A review of the work of Fabula and Hoyt also reveals that: (a) drag 
reduction effectiveness tends to improve with increased concentrations 
up to a certain concentration level, after which any further increase 
in concentration level will result in decreased effectiveness; (b) at 
concentration levels less than the “optimum" concentration for any 
given additive, the friction factor versus Reynolds number 

curves are roughly parallel to, but lower than, the friction factor 
curve for the pure solvent flowing in the turbulent flow regime; 

(c) at concentration levels above the optimum concentration the friction 
factor curves are roughly parallel to, but higher than, the extended 
laminar fluid flow curve of the pure Newtonian solvent. Additional 
effects which appear to be more related to the flow system than to the 
fluid are also commonly observed: 

(b) For a given solution which exhibits drag reduction 
characteristics the drag reduction effect is more pronounced as the 
diameter of the flow system is increased. 

(b) For a given solution and flow system, the drag reduc- 
tion effect becomes more pronounced as flow rates are increased. 
However, there appears to be a maximum tolerable shear rate associated 
with each flow system beyond which the drag reducing effectiveness 


appears to decline. 


The more recent work of Shin (25) has served to quantify some of 
these qualitative observations, particularly with respect to defining 
miemsignificant characteristics of a "dilute," “concentrated," "critical" 
and “optimum" solution. Shin's results suggest that the key to the 
question of whether the polymer solution is dilute or concentrated 
depends on the critical concentration which corresponds to the condition 


that each polymer molecule (having been assigned an effective spherical 








diameter which is of the order of magnitude of its r.m.s. end-to-end 
distance) is brought to a spherical packing of 74%. In other words, 
the random coils touch each other, but are otherwise unchanged. When 
the polymer concentration is less than the critical concentration, it 
is considered “dilute." It i1s%whenva solution 1s. "concentrated sunar 
it will display considerable gradient dependency of viscosity or non- 
Newtonian behavior. On the other hand, polymer solutions less than one 
third of the critical concentration appear to display either negligible 
or undetectable gradient dependency of viscosity. While there was no 
direct correlation achieved between the optimum concentration and 
critical concentration, the optimum concentration was observed to be 
generally one tenth to one hundredth of the critical concentration. 

For example, the critical concentration of polyisobutylene L-200 in 
cyclohexane was found to be 520 wppm; the optimum drag reduction effect- 
iveness over the limited shear rates investigated was found to be 

35 wppm. As was emphasized by Shin, it is not commonly recognized that 
the critical concentration of solutions of polymers having molecular 
weights in the millions is small. Indeed, 0.1% polymer solutions are 
usually "concentrated" as defined by Shin. 

In even more recent work, Rodriguez et al. (20) appear to have 
achieved a very useful correlation between drag reducing characteris- 
tics for turbulent flow in a pipe and measurable properties of the 
polymer solution and flow system. They have successfully correlated 
friction reduction behavior of a viscoelastic solution with a modified 


[beborah number" defined as: 


2 


Ne SAY t,/D [1-5] 


Db 

Fabula et al. (7) and more recently Hershey (8), Astarita (2) 
and others have suggested that fluid friction reduction occurs when the 
relaxation time of the polymer molecules in solution exceeds a character- 
istic flow time of the experiment. The relaxation time of the polymer 
solution represents the relative amounts of viscous and elastic response. 
Characteristic flow time has been taken as the reciprocal of the shear 
rate at the wall. The problem has been that of identifying specific 


physical quantities susceptible to direct measurement, which would be 
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characteristic of fluid relaxation time. Park (13) has shown that- 
intrinsic viscosity will not yield a favorable correlation when evalu- 
ating concentrated polymeric solutions. Yet the basis for determining 
a measure of relaxation time (t,) in the Rodriguez correlation would 
appear to be intrinsic viscosity (though fluid velocity and flew system 
diameter are, of course, reflected in the modified Deborah number). 
There would appear to be some disparity between the results of Park and 
Rodriguez. It would be interesting to investigate the drag reducing 
capabilities of Carbopol using the Rodriguez correlation since this 
additive would appear to be ideally Suited to promote drag reduction, 


though none has been observed in Carbopol solutions to date! 


ee a em ee ee ee ee 


The original scope of this investigation envisioned an examina- 
tion of the drag reducing effects of a number of polymers in several of 
the more commonly used military fuels. It was further recognized that 
the effect of these additives would have on the quality of the fuel 
would be of paramount interest. However, it would first be necessary to 
design, assemble and calibrate a flow system which would be suitable 
for investigation of the drag reduction phenomena. Because of the time 
involved in this phase of the work, the original scope was limited to 
a brief examination of a high molecular weight additive in military jet 
fuel at a concentration of 50 weight parts per million. 

Time also placed important limitations on the type of flow system 
which could be assembled and used for drag reduction experiments. 
Although a closed loop, recirculating type of flow system, patterned 
after those of Dodge (4) or Melton and Malone (12) was initially 
considered, it became clear that the time required to obtain aT assemble 
the equipment needed for such an elaborate system would be excessive. 

It might also be pertinent to note that systems which provide for a 
recirculation of the polymeric solution will not permit as detailed an 
examination of such important factors as mechanical degradation as is 
permitted in a single pass flow system. This is particularly true if 
recirculation is accomplished by a mechanical device generating high 


localized shear rates. For these reasons, it appeared that a single 
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pass system patterned after a variable pressure capillary tube visco- 
meter (such as the one used by Dodge for viscometric measurements) 
would be the most practical system which could be assembled in a reason- 
able time period. Hypodermic tubes for use as test sections were readily 
available as was a supply of bottled nitrogen gas for use as the 
driving mechanism. Since it was desirable to design a flow system 
capable of providing a means of observing the flow behavior in the well 
developed turbulent regime, a computerized model of a capillary tube 
viscometer was developed through which the design parameters affecting 
fluid flow over a range of Reynolds numbers up to 10° could be examined. 
The maximum driving pressure, of course, could not exceed the pressure 
provided by the bottled nitrogen. Further, the largest fluid reservoir 
avallable capable of withstanding the pressure from the bottled nitro- 
gen appeared to be an empty nitrogen cylinder which had an estimated 


capacity of between 10 and 15 gallons. 
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DESIGN, CONSTRUCTION... CALIBRATION. AND OPERATION OF FLOW SyeiEn 


PRELIMINARY DESIGN CONSIDERATIONS 


Previous experimental work involving observations on drag reduc- 
tion have generally been performed on three distinctly different types 
of flow systems: 

(a) The immersed rotating disk or cylinder. 

(b) The closed loop, recirculating flow system. 

(c) The open loop, single pass flow system through small 
| capillary tubes. 


Generally, these involved a type of rotational shearing boundary 
or cylindrical tube flow system, each of which has its own distinct 
advantages and disadvantages. Thts study was undertaken with a specific 
application in mind--the feasibility of reducing the time required to 
perform certain critical refueling operations encountered in the logis- 
tics support of tactical military maneuvers. A characteristic of these 
operations is the "one shot" transmission of fluid through relatively 
short fuel lines. A modification of the latter system was selected for 
the experimental work in this study in order to achieve a degree of 
similarity to the characteristics of tactical military logistical 
systems. 

The system consisted of a set of stainless steel hypodermic tubes 
which could be connected to a high pressure reservoir into which the 
fluid to be examined was loaded. The fluid was forced out of the reser- 
voir by a siphoning technique which was achieved by connecting the 
reservoir to a source of bottled nitrogen (see Figure 2-1). 

An energy balance on a capillary tube flow system reduces to the 
following form, considering steady state, isothermal flow and neglecting 


the effects of hydrostatic pressure: 
2 Z 2 _ 
(Poe) /e + (K/2g )(V, ae ) + ea L/g D a [2-1] 


where points 1 and 2 refer to the entrance and exit to the flow tube. 


a2 
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In capillary tube measurements, it is customary to measure the total 
pressure applied to the reservoir chamber (P,) and to consider the 
velocity in the reservoir (Vv, ) essentially as zero. Equation [2-1] 


then reduces to: 
KP = (pV"/2g.) Oe aly Ab) {2-2 | 


The applied pressure then is used to overcome viscous resistance to flow 
within the tube and to provide the kinetic energy of the stream as well 
as the entrance and exit losses. These latter forms of energy are 
usually accounted for in the variable K. For ordinary laminar Newtonian 
flow probably the most widely accepted value of K is 2.24 (4). Unfor- 
tunately, the appropriate value of the total correction to account for 
kaometic energy and entrance and exit effects is somewhat in doubt for 
non-Newtonian and viscoelastic fluids (24). One can of course, place 

a lower limit on the value, that of 1.0 for a kinetic energy correction 
corresponding to a perfectly flat velocity profile (thus, a correction of 
ev*/2g) Dodge (4) arbitrarily assigned a correction corresponding for 


the combined entrance, exit and kinetic energy effects of 1.5 pv" /g, 


corresponding to a value of K equal to 3. 

The ideal solution to this dilemma is to make the unknown value 
of K insignificant when compared to the 4fL/D term. In the laminar 
regime--for example at a Reynolds number of 800--f takes on a vlaue of 
G02. A length to diameter ratio of 1:1000 will result in a value of 
4UfL/D =80. Under these conditions, the widest conceivable error for K 
is probably less that 1%. If kinetic energy, exit and entrance effects 
are completely disregarded, the error would be less than 3%. On the 
other hand, in the turbulent regime with a drag reduction mechanism 
operative, the observed Fanning friction factors will frequently be less 
than 0.002. The same uncertainty in K as in the previous example will 
now result in an error of almost 8%. Further, it was found that an L/D 
ratio in excess of 200 could not be tolerated if the goal of observa- 
tions extending to a Reynolds number of fig was to be achieved for the 
full class of fuels originally envisioned. The uncertainty in K intro- 
duced by this new criterion could perpetuate an error of as much as 20%. 


Specially designed flow tubes with pressure ports located well downstream 


LS 
from the entrance to the tube obviated the need to consider these 
corrections. 

It might be pertinent to note that drag reduction experiments 
performed in a capillary tube rheometer at the University of Minnesota 
were reported in 1964 (19). The flow tubes were 1000 diameters in 
length in all cases. The report indicates that this “arbitrary selec- 
tion of length was based on the recognition that fluid characterization 
errors contributed by entrance, exit, etc., diminished with increasing 
length while errors contributed by shear degradation might increase with 
length." However, the highest Reynolds number achieved in this apparatus 
with drag reduction effects operative was approximately 25,000. 

A computerized model of a capillary tube flow system was developed 
through which Senne design concepts could be tested for feasibility 
and practicality. 

The computer model used for preliminary design was expanded as 
the final design of the actual flow system progressed, so that the effect 
of each design feature could be tested in advance of final construction. 
Therefore, computer program number 1 in Appendix A accounts for the 
design characteristics peculiar to the permanently installed portions of 
the flow system. This program can be used for two different purposes: 

(a) To provide the reservoir pressure which will be required 
to operate the flow system at a given Reynolds number. 

(b) To test the effect that variations of flow tube char- 
acteristics such as flow tube diameter, entrance length, exit length, and 


test section length will have on the performance of the flow system. 


The former feature is designed to assist the operator of the flow 
system by providing the minimum weight or volume of fluid required to 
meet the criterion that the flow time for the run last at least 100 
seconds, and by providing the pressure regulator setting required to 
achieve the flow regime (at a given Reynolds number) desired. The latter 
feature is provided to serve as a guide to the design and expected per- 
formance of new flow tubes to be used with this flow system. In either 
case, if the system limitations of available pressure or reservoir capac- 


ity are exceeded, the program will automatically signal such a condition 
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to the operator. A sample write out from this program is included in 


Table 2-1. 


DISCUSSION OF SUPPORT EQUIPMENT USED AND CALIBRATION TECHNIQUES. BME OWED 


re re mney er greg rere ee 


Fluid flow rates were determined by timing the flow of a pre- 
determined quantity of fluid with a Galco stop watch (Chemical Engineer- 
ing Department #110). The watch dial indicated time in increments of 
0.2 seconds but there was a sufficient spread between marked increments 
to allow interpolation to the nearest 0.1 second. A predetermined 
criterion was lmposed: that a run must take at least 100 seconds to 
insure that the precision of the time measurement was 1 part in 1000. 
The decision as to whether a volumetric or weight measurement would be 
taken using 100 and 500 ml graduated cylinders, a 2000 ml volumetric 
flask and a 4000 Erlenmeyer flask which was calibrated with the 2000 ml 
flask and marked at the 4000 ml level. When it was determined that the 
rate of flow would be too great to permit a time measurement of 4000 ml 
in at least 100 seconds, a weight measurement of fluid flow was taken 
using a Fairbanks Morse and Company 200 pound capacity scale, serial 
number G 403886, code 528. This scale contained marked increments of 
1/16 pound with a sufficient spread between the marked increments to 
permit interpolation to the nearest 1/32 pound. A fluid sample of at 
least ten pounds was used as a basis for determining the rate of flow 
when measurements were taken, using a ten gallon stainless steel con- 
tainer to collect the sample of fluid. Comparison of the volumetric 
and weight measurements used to determine flow rate, established by 
weighing a 4000 ml fluid sample of carefully measured density, indicated 
the percent difference to be 0.6%. This difference corresponded to the 
minimum weight discrimination believed achievable on the Fairbanks Morse 
Scale (1/32 pound). 

Fluid densities were measured with a Westphal balance obtained 
from the W. M. Welch Scientific Company. The weights were identified by 
serial number 4028A and the bulb was numbered 472 WHR. This balance was 
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MODEL OF EYPECTED FLOW PERFORMANCE 
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initially calibrated by using distilled water, the density of which has 
been established accurately (16). The balance was placed on a smooth 
aluminum plate in a semi-enclosed cabinet to minimize environmental dis- 
turbance. The calibration was checked periodically during the research 
and calibration phases to insure against calibration shift caused by 
accidental movement. Using this balance, densities were measured in 
increments of 0.0001 grams per cubic centimeter up to a maximum of 1.8999 
grams per cubic centimeter. 

Temperatures of both the environment and the fluids were taken 
using a Princo 15-30°C thermometer calibrated in increments of 0.1°C with 
a sufficient spread between marked intervals to permit interpolation 
meethe nearest 0.05°C. 

Three different pressure regulators connecting the nitrogen source 
cylinder with the fluid reservoir were used depending upon the range of 
pressures required to drive the fluid through the tubes. Regulated 
pressures of 200-1500 pounds per square inch were attained using a 
Hoke-Phoenix 0-3000 psi regulator, style 521 serial number N296C. 
Pressures of 75-175 psi were obtained through the use of a REGO 0-400 psi 
regulator, type XL. A Victor Equipment Company 0-60 psi regulator was 
used when pressures ranging from 10-60 psi were required. 

Viscosity measurements were taken with a calibrated Cannon- 
Fenske capillary tube viscometer number 25/J680. The supplied calibra- 
tion constant of the viscometer was 0.002225 centistokes per second at 
mOOCES and 0.002218 centistokes at 210°F. The calibration constant at 
70°F obtained through linear extrapolation of supplied values was deter- 
mined to be 0.002227 centistokes per second. 

A series of calibration runs to check this value was conducted 
using distilled water for which accurately measured data on viscosity 
and density as a function of temperature was readily available (16). A 
computer program (program number 2 in Appendix A) was then prepared 
which would accept values of measured time and temperature and the cor- 
responding known values of viscosity and density of the distilled water. 
From this information an average calibration constant was determined for 
the viscometer over the environmental temperature ranges encountered. 


ie results of these calibration. runs are included in Table 2-2. The 
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RESULTS ON THE CALIBRATION GF NO. 25/J680 
CAPILLARY TUBE VISCOMETER 


COMPUTED 
FLOW FLUID STANDARD STANDARD KINEMATIC CALIBRATION 
HIE TEMPS GENSITY  VISCOST) Ye Sees a yee ane 


Mees) (DEG.C) (GH/CC) (CP) (CS) (esc ecn 
431.4 21.6 0.9979 0.9669 0.9690 6.00224614 
431.6 Diem 0.9979 0.9669 0.9690 0,9022451 
432.8 21.8% 0,9978 0.9623 0.9644 0.0022283 
431.0 noo 77 0.9510 0.9532 0.0022116 
426.0 Pee Onoo7 7 0.9444 0.9466 0.0022220 
424.0 22.9 Oro 7G OTe 0.9400 (i 62 ae 


COMPUTED ARITHMETIC AVE. CALIBRATION CANSTANT 0.0022283 
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average calibration constant of 0.002228 was used throughout the exper- 
imental phase of this study as the basis for determining viscosity with 
the Cannon-Fenske viscometer. 

The pressure and differential pressure measuring instruments used 
in these experiments include a Heise 0-2000 psi bourdon pressure gage 
number H2939R, a Midwest Instrument 0-1000 psi differential bourdon 
pressure gage and a 50 inch glass manometer tube suitable for safe use 
at pressures corresponding to approximately two atmospheres. Both 
bourdon gages are marked to indicate increments of 5 psi over their 
entire range. The spread of the marked increments on the Heise gage is 
sufficient to permit reliable interpolation to within 1 psi while the 
design of the Midwest gage is such that a reliable een Ee 
Within 2 psi was the best accuracy believed attainable. 

The Heise and Midwest gages were initially calibrated against a 
precision dead weight pressure measuring instrument, Ruska Balance 
#11516 with piston #B3-212, in the range from 100-1000 psi. While the 
Heise gage proved to be in exact agreement (within the 1 psi limit 
attainable for this gage) over the entire range considered, the corres- 
ponding readings taken with the Midwest gage raised a serious question 
of repeatable reliability. In order to establish a statistical base from 
which to judge the repeatability and usefulness of this gage, a series 
of recalibration runs were conducted. The Heise gage and Midwest gage 
were connected in parallel to a source of bottled nitrogen through the 
Hoke-Phoenix pressure regulator and the pressure was systematically 
raised and lowered over a range from 0-700 psi. The Heise gage was 
used as a secondary standard against which the performance of the Midwest 
gage could be measured. The results of this calibration are indicated 
in Figure 2-2. A mechanical hysteresis of 2 psi over this pecotune 
range is readily apparent. 

While the statistical scatter was considerable in the lower pres- 
sure ranges, it was decided that the repeatability was sufficiently good 
above the 75 psi range to justify the continued use of this gage. How- 
ever an allowance providing for the observed statistical uncertainty has 
been included in all calculations derived from values measured with this 


gage. While it is recognized to be the measuring instrument with the 
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least desirable characteristics of all the equipment used, it was the 
only differential pressure gage available that was capable of measuring 
the high differential pressures expected to be encountered in the 
course of the experimental investigation. Further, time limitations 
would not permit the further delays which would be required to obtain, 


install and calibrate a different gage. 


CONSTRUCTION AND CALIBRATION OF THE VARIABLE PRESSURE FLOW SYSTEM 


mr a SSS Sa SSS 


A specially adapted high pressure nitrogen cylinder was employed 
as the fluid reservoir of the flow system. A safe working pressure of 
2200 pounds per square inch was determined by tests conducted in July 
of 1962. The capacity of the cylinder was determined to be 1.78 cubic 
feet (111.5 pounds of tap water measured at 17.9°C). The exit to the 
nitrogen cylinder was fitted with a specially designed cap capable of 
receiving nitrogen gas under high pressure and discharging the sample 
fluid to a flow tube. A detailed drawing of the cap is shown in 
Figure 2-3. 

One of the unique features of the flow system is the use of 
specially constructed tubes (detailed drawings of these tubes are shown 
in Figure 2-4). Hypodermic tubing with inside diameters varying from 
0.054 up to 0.106 inches was obtained from the C. A. Roberts Company. 
Five tubes were constructed each with test sections of approximately 100 
diameters, entrance lengths of 75 diameters and exit lengths of 25 
diameters. A small hole was drilled at either end of each test section 
to which 1/8 inch O. D. stainless steel tubing was connected through 
5/16 inch O. D. copper sleeves. The sleeves were reamed to snugly 
enclose the tubes, thereby providing support to the flow tubes as well 
as to the 1/8 inch tubing through which pressure measurements could be 
taken. The exit end of the flow tube discharged to a larger tube which 
was used to carry and deflect the efflux to a collection receptacle. 
The entrance end of the tube was fitted into a 3/8 inch Autoclave 
Engineering plug fitting which could then be assembled into an Auto- 
clave series 30, 90° angle valve. The entrance rim of the Autoclave 


- plug was rounded to meet the tubes thereby minimizing the disturbance 
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caused by an abrupt decrease in flow diameter at that particular point 
in the flow system. Hypodermic tubes with an outside diameter very 
nearly equal to the inside diameter of the flow tubes were used to 
deburr and polish the flow tubes after they had been constructed and 
during their experimental use. 

The length of the test section of the flow tube was measured to 
better than 1/64 inch using a precision machinist's ruler. However, 
since the internal diameter of the hydraulic tubing supplied by the manu- 
facturer was guaranteed to be accurate within limits of approximately 
+0.002 inches, it was necessary to determine this measurement more pre- 
cisely by conducting flow runs in the laminar flow regime. For this 
purpose tap water was loaded into the fluid reservoir and forced through 
the flow tubes at Reynolds numbers no greater than 2,000. A manometer 
containing colored carbon tetrachloride which had a measured density of 
1.557 grams per cubic centimeters was connected across the pressure 
ports of the flow tube. Since the difference in density of the mano- 
metric fluid and the tap water was approximately 0.56 grams per cubic 
centimeters, this device provided a very sensitive and accurate means 
of determining the pressure differential across the test section. 

A computer program was then prepared to compute the internal 
diameter of the flow tubes from the following measured values: (a) test 
section length in inches; (b) fluid density in grams per cubic centi- 
meter; (c) fluid viscosity as measured by the Cannon-Fenske viscometer 
im seconds; (d) volume of fluid flow in milliliters; (e) the time of 
Fluid flow in seconds; and (f) pressure in inches of manometric fluid. 


The basis for the computation is the well known Hagen-Poiseuille law: 
= a 
Das m(P. P,) ER /8uL [2-3] 
which appears in computer program #3 in Appendix A in the form of: 
4 
D = [(128uQL)/(rAPg ) J hea 


The results of these calibration runs are summarized in Table 2-3. Com- 
puter listings for each of the five flow tubes calibrated in this manner 


appear in Appendix B. 
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TABLE 2-3 
SUMMARY OF FLOW TUBE DIMENSIONS 


Manufacturers Test 
Flow Tube Supplied Calibrated Section Entrance Exit 
Gage Diameter Diameter Length Length Length 
Number (Inches ) (Inches) | (Inches) (Inches) (Inches) 
1s OeOou O20 Sco) S206 3.99 dea 
kS 00:71: 007 Se Tako So as ee 7s 
12 0.085 0.08548 SaO0 O35 Dace 
ae]! 0.094 0.09264 9.460 OS Doo 
10 0.106 0.10590 10.609 7.95 2.65 


In addition to providing a means by which the internal diameter 
of the flow tubes could be accurately measured, the results of these 
tests indicated that data of reasonable consistency could be obtained 
in the laminar flow regime using the measurement techniques and instru- 
ments employed. 

Performance characteristics of the tubes was next tested in the 
turbulent flow regime corresponding to a range of Reynolds numbers from 
10" to 10°. While it would have been desirable to measure these char- 
acteristics for each of the five tubes, time limitations permitted 
detailed examination of only two tubes, number 15 and 13. These calibra- 
tion experiments were conducted over the same turbulent flow regime 
using two fluids, military jet fuel (JP4) and tap water. The flow 
behavior of each of these fluids was thereby independently reproduced. 

There were three major objectives to be accomplished during 
these performance tests: 

(a) To determine the interior roughness characteristics of 
the flow tubes. 

(b) To determine the effects, if any, of the pressure ports 
placed at either end of the test section over the range of Reynolds 
stresses to which the fluids would be subjected during research experi- 


ments. 
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(c) To observe the performance of the Midwest differential 
pressure gage aS a measuring device under dynamic turbulent flow condi- 


TOT Se 


It might be well to recall, in connection with this latter 
objective, that the characteristics of the Midwest differential pres- 
sure gage were originally determined by a calibration technique which 
produced a statte pressure differential across the two inlet ports of 
the gage. Further, there was a clearly defined mechanical hysteresis 
associated with this gage which appeared when static pressures were 
being systematically increased and decreased. This observed static 
behavior raised serious questions of the gage's ability to yield reli- 
able pressure measurements under conditions of turbulent fluid flow 
which is characterized by violent pressure fluctuations. 

It was recognized, of course, that the effects of the pressure 
ports could only be judged if the interior roughness characteristics 
failed to conform with well established flow behavior data. The dilemma 
just described is not uncommon to scientific research: Experimental 
results cannot be achieved without measurement--yet to measure is to 
disburb. Further, whether a meaningful discrimination between the 
first two objectives could be achieved would largely depend upon the 
performance of the differential pressure gage. 

There were, however, some criteria by which the results of these 
tests could be judged: 

(a) The results of the calibration tests on the differential 
pressure gage suggested a decrease in reliability at the lower pressure 
ranges. Further, at pressure readings of less than 100 psi, the number 
of significant digits available was reduced by one. _ 

(b) Any abnormal effects on flow behavior resulting from the 
pressure ports could be expected to become more pronounced as the flow 
rates were increased. 

(c) The expected minimum friction factor at any particular 
flow rate would correspond to the friction factor established for 
"hydraulically smooth" pipes. The expected maximum would not be likely 


_to exceed the friction factor relationship established for drawn tubing. 
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For the 15 gage flow tube with a 0.053 inch I.D., the maximum expected 
relative roughness waS approximately 0.002. The corresponding value 


for the 13 gage tube was approximately 0.001. 


Using these criteria, the results of the turbulent flow tests 
as illustrated in Figure 2-5 suggest the following conclusions: 

(a) The reliability of pressure measurements taken with the 
Midwest gage is substantially reduced at differential pressure readings 
below 75 psi. 

(b) The relative roughness of the flow tubes, including any 
contributing effects of the pressure ports, is between 1/5 and 1/10 the 
relative roughness associated with drawn tubing. 

(c) No abnormal flow behavior can be specifically identified 
with the pressure ports over the range of shear rates examined during 


mime test. 


Hence, it was concluded that the equipment and flow system design 
were entirely satisfactory for experimental studies on drag reduction. 
A computer write out of information used to construct Figure 2-5 is 
contained in Appendix C. Computer program number 4, used to produce 


this information is contained in Appendix A. 


OPERATION OF FLOW SYSTEM 


Once a sample had been loaded into the fluid reservoir a choice 
of the flow tube to be used and the range of the flow regime to be 
investigated was required in order to determine appropriate measurement 
equipment. Generally, the Midwest differential pressure gage was 
employed to measure the pressure drop across the test section at high 
Reynolds numbers and the flow rate was measured on the Fairbanks Morse 
Scale by timing a predetermined wetght of fluid. The mercury manometer 
was employed to measure the test section pressure differential at the 
lower Reynolds numbers and the flow of a predetermined volwne of fluid 
(up to a maximum of 4000 milliliters) was timed to measure flow rate. 

The computer flow model commented upon earlier in this Chapter 
was used as a guide for these determinations. By referring to Table 2-l, 


“for example, if the flow of pure tap water through the 15 gage tube at a 
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Reynolds number in the vicinity of 60,000 was to be observed, the pres- 
sure regulator from the nitrogen cylinder was set at 800 pounds per 
Square inch and the program indicated that the flow of at least 15 
pounds (rounded to the nearest 5 pounds) of the fluid should be timed 
in order to insure that the timed run would consume at least 100 sec- 
onds. Therefore, a series of runs could be preplanned to maximize the 
efficient use of fluid available for investigation. 

Once the flow regime to be investigated had been determined and 
the appropriate measurement equipment assembled, it was next necessary 
to purge the tubing leading from the pressure ports of the flow tubes 
to the pressure gages. This was accomplished by loosening the fittings 
at the entrance to the gages, plugging the exit end of the flow tube 
and permitting the fluid to flow through the pressure gage line until 
a free flow, clear of gas bubbles was observed. Next, a 500 milliliter 
sample of the fluid was withdrawn from the reservoir from which the 
basic properties of temperature, density and viscosity were determined. 
An average measurement of viscosity was obtained on the sample by 
timing three to five consecutive trial runs in the Cannon-Fenske capil- 
lary tube viscometer. Since there was not a means provided to maintain 
the reservoir fluid at a constant temperature, the fluid properties of 
the reservoir mixture were measured every 20 to 30 minutes during the 
course of experimentation. Variations in these properties for a given 
reservoir load were generally found to be negligible since the maximum 
ieee required to empty the cylinder rarely exceeded 45 minutes once 
experiments were begun. The cylinder used with the Westphal balance and 
the capillary tube viscometer were each cleaned and dried immediately 
after use. 

Prior to a test run, the weights on the Fairbanks Morse Scale 
were adjusted to permit two to four pounds of the test fluid to flow 
before the stop watch was started and measurements were begun. This pro- 
cedure allowed a sufficient amount of time for the Autoclave 90° valve 
to be fully opened and permitted the fluid flow pattern to develop and 
stabilize before measurements were taken. When the arm of the Fairbanks 
Morse Scale returned to equilibrium, the stop watch was started and the 


predetermined weights were added to the balance arm. Measurements of 
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gage pressure at the upstream end of the test section and differential 
pressure across the test section as well as the weight of fluid to be 
timed were recorded on a data sheet. When enough fluid entered the 
collection container on the scale to cause the arm to return to equil- 
ibrium, the stop watch was stopped and the flow valve was closed. The 
temperature of the fluid was then taken and recorded together with the 


elapsed time for the run. 


ANALYSIS OF FLOW DATA 


One basis for determining the effectiveness of a polymer would 
be to relate the friction factor observed for a fluid with such an 
additive to the friction factor observed for a purely viscous flutd 
flowing at the same Reynolds number. The factors affecting the Reynolds 
number--ie., density, average fluid velocity, tube diameter and fluid 
viscosity--could all be accurately determined to three or more signifi- 
cant figures. However, the friction factor varies as a direct function 
of pressure which, when measured with the Midwest differential pressure 
gage, could be determined to two significant figures at best. Informa- 
tion indicating the scatter in pressure readings expected over the wide 
range of pressures investigated was available from the calibration runs 
performed on the gage (see Figure 2-2). This scatter includes an allow- 
ance for both the observed mechanical hysteresis as well as the random 
scatter associated with other limitations of the instrument and observer. 
This information was used to compute error bounds for each value of 
friction factor, assuming that the measurements taken with the differen- 
tial pressure gage represented the least reliable data and would 
contribute the most significant error. 

At the conclusion of a series of runs, data cards were prepared 
using the information taken during each run. These were used with 


Computer Program number 4 in Appendix A to analyze the test data. 


PREPARATION OF POLYMER SOLUTIONS 


The pioneering research undertaken by Fabula and Hoyt (6) at the 


Naval Ordnance Test Station during the late 1950's and early 1960's 
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was instrumental in identifying a relatively large number of water 
soluble polymers capable of inducing drag reduction effects in extremely 
dilute aqueous solutions. However, an equally extensive program 
directed toward the identification of additives capable of producing 
drag reduction effects in dilute, non-polar hydrocarbon solutions has 
never been undertaken. In fact, prior to Shin's investigation (23), 
Toms original research (25) appears to be unique, not only because it 
represents the first written account of drag reduction associated with 
dilute Newtonian solutions, but also because his experiments involved a 
non-aquecus solvent. Shaver and Merrill worked with concentrated solu- 
tions of polyisobutylene in cyclohexane (22). During World War II, an 
apparent reduction of flow resistance was noted when gasoline was 
thickened with napalm for use in flame throwers (1). However, in 
contrast to the dilute solutions studied by Toms and Fabula,. the final 
solutions involved in these isolated examples had markedly different 
fluid characteristics than the criginal solvents. The fact that poly- 
isobutylene had been shown to be soluble in a pure non-polar hydro- 
carbon solvent and the fact that it apparently displayed anomalous 
behavior in concentrated solutions clearly suggested that further study 
of its behavior in dilute solutions might be productive. Additionally, 
it seemed appropriate to test the more common water soluble polymers for 
solubility in jet fuel. For this reason, samples of several polymers 
known to exhibit viscoelastic properties were initially requested and 
obtained from various industrial sources. While most of these polymers 
were known to be soluble in water and sea water, the solubility of each 
was tested in the jet fuel for possible further examination if time would 
permit. A summary of the polymer samples and the sources from which they 
were obtained is contained in Table 2-4. 

Several 50 weight part per million solutions in jet fuel were pre- 
pared, each containing one of these polymers, and all were observed for 
several days. Each solution was prepared by carefully weighing a small 
sample of the polymer (between 34 and 38 milligrams) on a precalibrated 
precision Metler analytic balance and then adding the volume of jet fuel 
of carefully measured density that would. produce a concentration of 50 


weight parts per million. These samples were pleced in a one liter 
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TABLE 2-4 
LIST OF PCLYMERS OBTAINED AND INDUSTRIAte Sel ee 


Company Product 
Archer Daniels Midland Polysacchoride 
Minneapolis, Minn. ADM 7079 
Enjay Compeny Polyisobutylene 
New York, N.Y. L80, L140 and L200 
Union Carbide Corpcration Folvex WSR P301, 

New York, N.Y. Wek E205. 
N3000 and N750 

Hercules Powder Company Reten 210, Al and AS 

Wilmington, Del. 

So ee lecmee... Ine. Jaguar J2S1 

New York, N.Y. Polynal 27 


bottle and were periodically tested for solubility by visual observation. 
None of the water soluble polymers appeared to be soluble in jet fuel when 
mixed by natural diffusion techniques at room temperature. Although mix- 
ing by diffusion can be a very slow and difficult process for many poly- 
mers, it did not appear that more vigorous artificial methods should be 
employed in view of the well known susceptibility of most polymers to 
mechanical, thermal and/or chemical degradation. On the other hand, the 
three samples of polyisobutylene prepared in one liter bottles at con- 
centrations of 50 weight parts per million appeared to be completely dis- 
solved through natural diffusion when left standing at room temperature 
for approximately 36 hours. The preparation of more concentrated samples 
was found to be somewhat slower. However, even a 6% weight solution of 
the highest molecular weight polymer prepared (VISTANEX L-200) was 
completely dissolved as judged by visual inspection of such factors as 
homogeneity, clarity, etc., in less than two weeks. The speed with which 
the mixing occurred appears to be more a function of the total quantity 


of polymer present in a sample than the weight concentration of the 
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sample being prepared. That is, a 3% solution prepared in 500 ml size 
samples dissolved more rapidly than a 3% solution prepared in 1000 ml 
size batches. 

The concentrated polymer solutions were developed to facilitate 
the preparation of the 50 weight part per million solutions to be 
loaded into the fluid reservoir for flow experiments. It was found that 
these solutions once the polymer was fully dissolved, mixed very readily 
with the pure solvent to form the desired dilute concentration. Thus, by 
pouring a solution containing 2.267 grams of the polymer into a 20 gallon 
stainless steel tank and then adding enough jet fuel to bring the total 
weight to 100 pounds, a 50 wppm sample could be prepared in a very few 
Minutes. In an effort to insure that all of the polymer originally 
prepared in the one liter bottle was actually included in the 100 pound 
solution, the bottle was flushed several times with the pure jet fuel 
which was used to form the final solution. 

Since the preparation of large quantities of even dilute concen- 
trations can become a problem, the procedure adopted in this study is 
certainly recommended. Further, certain physical properties which are 
impossible to describe quantitatively such as homogeneity, purity, etc. 
become more difficult to control as the size of the sample increases. 

Once the desired dilute solution was formed, it was gently poured 
into the fluid reservoir at a rate of about 1/2 gallon per minute to 
minimize any effects of mechanical degradation or depolymerization prior 


to actual flow through the test section. 
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EXPERIMENTAL RESULTS AND DISCUSSION 


All drag reduction experiments were accomplished by observing the 
flow behavior of solutions containing 50 weight parts per million of the 
high molecular weight Vistanex L-200 in military jet fuel. The two flow 
tubes with the smallest diameters were used: number 13 with an 0.0716 
inch I. D. and number 15 with an 0.0533 inch I. D. The flow system was 
primarily operated in the turbulent regime corresponding to a range of 
Reynolds numbers from approximately 10,000 to 100,000, although a few 
observations were also made on laminar flow. Values of observed friction 
factor and Reynolds number were calculated for each run using the same 
computer program (number 4 in Appendix A) as was used to analyze flow 
data pertaining to the pure solvents. 

The results of the turbulent flow experiments are illustrated in 
Figure 3-1. (The data is presented in Appendix D.) The friction factor 
versus Reynolds number relationship found to exist for purely viscous 
fluids flowing in the turbulent regime, as well as an extension of 
laminar flow relationships have been included in order to facilitate 
visual comparison. Error bands have been included on all data calculated 
from pressure measurements taken with the Midwest D.P.I. gage in order to 
illustrate the statistical uncertainty observed at the time the gage was 
calibrated. Flow lines have been drawn which place greatest reliability 
on the data points calculated from pressure measurements taken with the 
mercury manometer. These data points do not have error bands. 

In examining Figure 3-1 it should be noted that: 

(a) At any given Reynolds nimben, the drag reduction ¢efteanr 
appears to be more pronounced as the internal diameter of the flow tube 
increases. 

(b) Over the range of Reynolds stresses covered by these 
experiments, the drag reduction effect appears to first increase then 


decrease with increasing Reynolds numbers. 


The former occurrence is well known and has been reported since 
Toms original findings as the "diameter effect." Its importance, how- 


ever, is often overlooked. It is because of this effect that dynamic 
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Fig. 3-1. Log-Log Plot of Observed Friction Factor vs. Reynolds 
Number for the Flow of 50 wppm PIB in Jet Fuel 
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Similarity in the turbulent flow regime is. lost. The Reynolds Number 
Similarity Law for turbulent flow of fluids no longer applies to this 
unique class of fluids. On the other hand, this anomalous behavior 
does not occur in the laminar region of fluid flow. A friction factor 
equal to or slightly higher than the relationship which correlates the 


Fanning friction factor with the Reynolds number in the laminar regime: 
f = 16/N,. [3-1] 


was found to apply to this ¢lass of fluids. (This data is not illus- 
trated on Figure 3-1 which has been scaled so as to permit graphical 
amplification of the more important region from Npg?l0 However, the 
Statement can be verified by reference to the complete data included in 
Appendix D.) 

The latter occurrence is somewhat more unique. In general, it 
has been found that the drag reduction effect becomes more pronounced 
at increasing flow rates. Many experimenters (6), (18), (19), (23), 
have observed that there is a critical flow rate beyond which the drag 
reducing effectiveness of an additive appears to deteriorate. Precisely 
what determines this critical flow rate has not yet been determined. 
The critical flow rate appears to be very dependent upon polymer con- 
centration, temperature, the particular polymer used etc., although 
these observations are not confirmed in this study since only one poly- 
mer at one concentration was examined at ambient temperatures. This 
decrease in friction reducing effectiveness was found to occur at 
around a Reynolds number of 50,000. It is believed to be indicative 
of severe mechanical degradation. 

At this point in the experimental program, it was believed that 
greater insight into the drag reduction phenomena could be achieved by a 
study of this apparent degradation than would be obtained by continuing 
on the plan originally established for the study. This departure 
appeared appropriate for a number of reasons: 

(a) The recent work of Shin, (23), Hershey, (8), Rodriguez 
et al. (20) had served to confirm the drag reducing behavior of dilute 
concentrations of polymers in hydrocarbon solvents. The effects of 


varying concentrations and flow tube diameters were well documented. 
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(b) Any practical application of the phenomenon in the 
field of military logistics as envisioned, must consider the limita- 


tions imposed by mechanical degradation. 


This is not to suggest that these limitations cannot be advanta- 
geously employed. Since the particular additive employed in this study 
consists of a hydrocarbon polymer of paraffinic structure (26), some 
form of artificial degradation after fuel transmission has been 
accomplished might well be preferred as a means of overcoming fuel 
quality deterioration. | | 

Two methods of artificially imposing mechanical degradation were 
employed: (a) ©The first method involved the reuse of a solution which 
had already been subjected to experimentation. The degradation imposed 
in this manner might well be representative of the effects to be expected 
as a result of long distance transmission; (b) The second method 
involved a modification of the flow system. The modification intro- 
duced a means of imposing a very high localized shear point prior to 
the flow tube. It will be recalled that care was exercised in the 
original design of the system to eliminate the occurrence of excessive 
fluid shearing prior to the point at which the fluid entered the flow 
tube. Even the entrance to the flow tube was tapered to minimize the 
pressure loss at that location. Indeed, when the 15 gage flow tube was 
used to observe fluid flowing at a Reynolds number of 100,000 in the 
tube, the maximum Reynolds number achieved at any other location in the 
system was approximately 15,000. 

The modification for the second method of imposing mechanical 
degradation was accomplished by inserting a needle valve into the system 
just upstream of the flow tube. A very high pressure drop could be 
imposed across this needle valve to produce an extremely high localized 
Reynolds number. The degradation imposed in this fashion would tend to 
simulate the mechanical devices often found in a complete closed loop 
flow system which produce high localized Reynolds stresses such as 
centrifugal pumps, etc. This system is referred to as the "modified flow 


system." 
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The results of these experiments are illustrated in Figure 3-2. 
The mechanical degradation imposed by both methods employed is clearly 
apparent. The more severe form would appear to accompany the use of 
the needle valve in the flow system. The polymer appears to be so 
degraded that the fluid no longer exhibits any substantial drag reduc- 
ing characteristics. Even the solution used for the first time through 
this arrangement (labeled as "first pass through modified flow system") 
appears to be substantially more degraded than the solution being 
recycled through the flow system for the second time. These results sug- 
gest that flow system design can have a very important effect on the 
performance of a drag reducing additive. That is, 1f maximum performance 
from a drag reducing additive is to be achieved, mechanical devices 
which produce high localized Reynolds stresses should be avoided. Addi- 
tives should be injected down stream from such a device. 

The mechanical degradation associated with either method appears 
to have permanently destroyed some of the drag reducing potential of the 
additive. That is, the friction factors observed for the solution which 
was passed through the system for the fourth time were consistently 
higher than those observed during even the third pass. 

Further, there would appear to be a striking similarity between 
the effects on friction drag resulting from the use of a degraded solu- 
tion and what has been termed “the diameter effect.' This suggests that 
the diameter effect is in fact more properly identified as a degradation 
effect. That is, the friction reducing capability of a solution is not 
necessarily decreased by virtue of fluid flow performed in smaller tubes. 
Rather, the mechanical degradation has been increased. This is not 
unexpected since the mechanical shear applied to a fluid flowing at a 
given Reynolds number is substantially increased as tube diameter is 
decreased. Nonetheless, the friction reducing characteristics of the 
additive, even in a degraded state, continue to be pronounced when com- 
pared to the flow behavior of the pure solvent. 

Finally, note the trend of the flow behavior patterns developing 
at Reynolds numbers between 10,000 and 50,000 as illustrated in Fig- 
ure 3-2. The friction factor curve produced from experiments uSing a 


fresh 50 wppm polyisobutylene solution appears to approach an extension 
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Number for the Flow of 50 wppm PIB in Jet Fuel 
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of the laminar flow curve. On the other hand, the most highly degraded 
solutions almost coincide with the turbulent friction factor curve. 

These characteristics were also observed in a qualitative sense 
during the flow experiments. During the flow of tap water, pure jet 
fuel or jet fuel with highly degraded polyisobutylene at well developed 
turbulent flow rates, a very marked pressure fluctuation was observed 
which was indicated by a pronounced "flutter" of the Heise gage pointer 
and manometric fluid. This flutter could be interpreted as variations 
in the absolute pressure as measured at the wall of the tube correspond- 
ing to the turbulent fluctuations known to exist. On the other hand, 
as the less degraded solutions were used there was a noticeable decrease 
in the magnitude of the flutter associated with each pressure indicator. 
When the freshly prepared solution was used, no noticeable flutter 
could be visually detected. It is recognized that the measuring instru- 
ments employed did not possess the sensitivity required to reduce these 
qualitative observations to quantitative data. Nonetheless, the 
general pattern just described could be visually detected even if accu= 
rate measurements could not be taken and the qualitative significance of 
these observations should not be ignored. 

A practical explanation of this observed behavior might be found 
in different modes of flow of the solvent and solute. Although the flow 
of the solvent is known to be turbulent at the Reynolds numbers investi- 
gated, little is known concerning the mode of flow of the dissolved 
solute. It would appear that the combined effect of solvent and solute 
flowing at a given Reynolds number produces a unique flow regime, the 
nature of which depends upon the relative (or superimposed) contribution 
of each. The flow characteristics of the fresher, less degraded solu- 
tion, least resemble the turbulent characteristics of the pure solvent. 
The inference is that the pure solute contributes a stabilizing influence 
more typical of laminar flow, producing a "“pseudo-laminar" regime when 
the effect of the solute is found to dominate. The maximum stabilizing 
contribution of the solute could be expected to occur at what has been 
called the “optimum" concentration. At concentrations less than 
optimum, the turbulent contribution of the solvent begins to dominate. 


At concentrations greater than optimum the contribution of the solute 
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continues to dominate. However, because of the thickening action of 
the solute, the viscous drag could be expected to be somewhat greater 
than that obtained at an optimum concentration. 

Mechanical degradation occurring during a flow experiment can 
be expected to alter the relative influence of the solute by effectively 
changing the original concentration. Thus, a highly degraded solution 
originally containing 50 wppm of solute could be expected to behave like 
a much weaker undegraded solution. Although this explanation does not 
attempt to describe the mechanism responsible for drag reduction 
associated with the solute, it suggests that any attempt to correlate 
properties of solutions which exhibit drag reduction characteristics 
must consider the very important influence of mechanical degradation. 

Unfortunately, until such a correlation is found, there is no 
mathematical device through which the performance of a given flow system 
can be predetermined since the dynamic similarity laws of purely 
viscous fluids do not appear to apply to this unique class of fluids. 
However, it can be shown that the rate of fluid flow in a given flow 
system is proportional to the reciprocal of the square root of the 


iarecion factor. That is, 
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Or. Save (f/f? a= oR 


where: Q represents flow rate; 
f represents friction factor; 
v refers to a purely viscous fluid; and, 
ve refers to a viscoelastic fluid. 

Although the precise value to be assigned to fe is not known, 
this relationship is useful in assessing the potential benefits to be 
derived from the drag reduction effect during tactical military replen- 
ishment. For example, a 65% reduction in friction factor from a value 
of 0.00517 to 0.00181 was the highest achieved in this study. The capac- 
ity of the flow system to transmit fluids was thereby increased by 
approximately 69%. Reductions in friction factor of 80% are not 
uncommon (6, 23). If this reduction could be achieved at a Reynolds 
number of Oe the capacity of the system would be increased by approxi- 


mately 224%. In other words, the time required to transmit fuel under 
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these circumstances could be better than halved. Even the solution 
passed through the system for the fourth time during this study 
exhibited a 28% decrease in frictional drag at a moderate Reynolds num- 
ber of 56,000. The flow rate was thereby increased some 17.5%. 

Time did not permit the evaluation of quality deterioration 
introduced by the addition of 50 wppm of polyisobutylene in military 
jet fuel. It was always recognized that the application of this con- 
cept to military logistics could not proceed until this very important 
phase had been accomplished. However, the benefits of achieving sub- 
stantial drag reduction with relatively small amounts of additives are 


believed to be sufficiently attractive to warrant this further effort. 





CHAPTER IV 


CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 


The results of the experimental study on the flow behavior of 
polyisobutylene dissolved in jet fuel indicate the following conclusions 


and recommendations for future work. 


1. Performance of Filow System The) variable pressure, Singleup a 
flow system designed and constructed for use in this study was found to be 
particularly well suited for experimental studies on viscous drag reduc- 
tion. ‘Flow behavior up to a Reynolds number of 10° could be observed using 
the apparatus as constructed for this study, yet fluid consumption could 
be minimized by the proper selection of flow tubes. The single pass type 
test system facilitated the study of shear degradation since the fluid 
makes only one pass through the test section so that exposure of the 
fluid to shear stress could be closely controlled. Also, this type of 
system is quite similar to that used in tactical military refueling 
operations. 

The apparatus was found to perform well in both the laminar and 
turbulent flow regimes when tested with liquids of known fluid properties. 
No abnormal flow behavior was identified with the pressure ports of the 
smaller tubes (numbers 15 and 13) over the range of shear rates examined 
during the calibration tests. The larger tubes were not used for drag 
reduction flow tests. However, from an examination of the data taken to 
determine the internal diameter of the larger tubes, it would appear 
that the pressure ports might well be influencing the flow behavior in 
these tubes. Since no drag reduction tests were made using the larger 
tubes, an analysis of this effect was not performed. 

The Midwest differential pressure gage used to measure high 
test section differential pressures did not provide the precision 
or range required in this study. Its continued use in future work is 
not recommended. There are other differential pressure gages available 
which will provide better precision for a given range of pressure. For 


future work, several precise differential pressure gages, each selected 


uy 
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tO cover a part of the total pressure range anticipated ye. oemoon 
and 100% of maximum range, should be used. 
It is recommended that a means of controlling the temperature 
fluctuations of the fluid in the reservoir be provided. These fluctuations, 
caused by ambient temperature changes in the laboratory, made it necessary 


to use average values of density and viscosity. 


2. Performance of Vistanex L-200 as a Drag Reducing Addiiryemeiae 
results of this study indicate that a marked reduction in viscous drag 
can be achieved by the addition of minute quantities of high molecular 
weight polyisobutylene in a hydrocarbon solvent. Frictional drag de- 
creases of as much as 65% were noted when turbulent flow experiments were 
performed with a 50 wppm concentration of Vistanex L-200 dissolved in. 
military jet fuel. The friction reduction was observed to increase with 
increasing flow tube diameter. Further, friction reduction increased 
with increasing Reynolds numbers up to a certain "critical" Reynolds number, 
after which friction reduction decreased with increasing Reynolds numbers. 
By way of contrast, this anomalous behavior was not observed in the 
laminar region of fluid flow where frictional drag was slightly higher 


than that predicted by the Hagen-Poiseuille relationship. 


3- sensitivity of Vistanex L-200 to Mechanical Degradation The 
solutions of polymer tested were shown to be extremely sensitive to 
mechanical degradation. Further, repeated testing of a fluid containing 

a fixed amount of polymer indicated that the drag reducing characteristics 
of such a solution could be irreversibly destroyed. The degradation was 
found to occur when the solution was subjected to a region of highly 
localized shear (a needle valve) as well as to the continuous wall shear 

in the experimental tubes employed. 

Any practical application of the drag reduction phenomenon must 
consider the limitations imposed by degradation. Therefore, if operational 
use is contemplated in a full scale flow system, it is recommended that 
mechanical devices producing high Reynolds stresses be avoided or, alterna- 


tively, injection of the additive be accomplished downstream from such a 


device. 
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4. Relationship Between Diameter Effect and Degradation Although 
the results of this study are not conclusive, there appeayouuosboaas 
striking similarity between what has been termed "the diameter effect" 
and the effects of mechanical degradation. It is therefore recommended 
that additional work be undertaken to investigate the relationship 


between the flow tube diameter effect and the mechanical degradation 


effect of drag reduction. 


>. Pseudo-Laminar Fllow Characteristics Of Fresh Solutions) Were 
were strong indications in the experimental work, i.e. the pulsation of 
the pressure measuring devices, that a suppression of turbulent fluctu- 
ations might well be one of the major characteristics associated with the 
drag reducing influence of an undegraded polymer solution. However, it 
is recognized that the pressure measuring equipment employed in this 
study could not permit a quantitative verification of this observation. 
It is therefore recommended that additional work be undertaken to verify 


Mitese Qualitative observations with pressure sensors with high frequency 


response, e.g., strain gages. 
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NOMENCLATURE 


Definition 


derivative of .. 

flow tube internal diameter 
Fanning friction factor 
gravitational constant 


correction to pressure head to account for exit, 
entrance and kinetic energy losses 


PGW aEUDe sre sir ssccrronml cme 

parameter in power law model 

index of non-Newtonian behavior in power law model 
Deborah number 

Reynolds number 

pressure 

volumetric flow rate 

radius of flow tube 

average velocity 

instantaneous velocity 


cartesian coordinates 


Greek Symbols 
Newtonian viscosity 
Oe HAS ss 
fluid density 
shear stress 
first mode relaxation time of a viscoelastic solution 

Subscripts 

evaluated at point 0,1,2... 
purely viscous fluid 
viscoelastic fluid 
evaluated at the flow tube wall 


Garresian Coordinates 





10. 


ql. 


ere 


13. 


14. 
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COMPUTER PROGRAM NUMBER 1 
Computer Model of Newtonian Fluid Flow in Single Pass Variable 


mm i yp || ee eet 


Pressure Flow System. 


This program was developed to aid in the operation of the single 
pass variable pressure flow system used in this study. It can also be 
used to test the effect that variations in flow tube dimensions such as 
flow tube diameter, test section length, etc., will have on the performance 
of the flow system. A sample of the output to be expected from this 
program is found in Table 2-1. 

Input data is entered in seven fields each containing 10 columns 
team F10.0 format. The Reynolds number (RN) for which performance 
characteristics are desired is entered in the first field. If the 
number 200 is entered in this field, the program will compute and write 
out flow characteristics for fluid flow at various Reynolds numbers in 
the laminar and turbulent flow regimes up to a maximum Reynolds number 
of 140,000. 

The flow tube diameter (D) in inches is entered in the second 
field. Fluid viscosity (VISC) in centipoise is entered in the third 
field. Fluid density (SG) in grams per cubic centimeter is entered in 
the fourth field. The flow tube exit length (CLD), entrance length (CLU) 
and test section length (L), all in inches are entered in the fifth, sixth 
and seventh fields. If any one of these fields is left blank, the program 
will assume an entrance length of 75 diameters, exit length of 25 
diameters, or test section length of 100 diameters for the blank field. 

The program will assist in operating the flow system at a given 
Reynolds number by providing the minimum weight or volume of fluid 
required to meet the criterion that the flow time for a run last at 
HeaamelOO seconds and by providing the pressure regulator setting 
required to achieve the flow regime desired. If the system limitations 
of available pressure or reservoir capacity are exceeded, the program 
Will automatically signal this condition to the operator. A blank 
card at the end of the input data will exit the program from the computer. 


Sample input data are located at the end of the program. 
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COMPUTER PROGRAM NUMBER 2 
Computer Program to Compute a Calibration Constant for a Capillary 


me em ee ee 


imMibe Viscometer. 


This program was developed to compute an average calibration 
constant for a Cannon-Fenske Capillary Tube Viscometer from data taken 
during calibration runs. A sample of output to be expected from this 
program is found in Table 2-2. 

Input data 1s entered in four fields each containing 10 columns 
in an F10.0 format. The efflux time (TIME) in seconds of the fluid 
flow in the capillary tube is entered in the first field. The fluid 
temperature (TEMP) in degrees Celcius is entered in the second field. 
The fluid density (RHO) in grams per cubic centimeter and fluid vis- 
cosity (VIS) in centipoise corresponding to the measured temperature 
are entered in the third and fourth fields. The program will compute 
a calibration constant for each calibration run and the arithmetic 
average of the series of runs. A blank card at the end of the input 


data deck will exit the program from the computer. 
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COMPUTER PROGRAM NUMBER 3 


Computer Program to Calculate the Internal Diameter of Flow Tubes. 


ee ae creas er ree 


This program was developed to compute an average internal diameter 
for the flow tubes to be used with the single pass variable pressure flow 
system used in this study. A sample of the output to be expected from 
this program is illustrated in Appendix B 

Each series of data cards containing information taken during 
flow experiments pertaining to a tube to be calibrated is preceded by 
a run control card. The run control card contains three data fields 
Beem of which contains 10 columns in an F10.0 format. The identifying 
number (such as the gage number) of the tube is entered in the first 
field (TUBE). The diameter (SPEC) of the tube in inches as supplied 
by the manufacturer is entered in the second field. The temperature 
(TEMP) of the calibration fluid is entered in the third field in degrees 
Celcius. 

The data cards following the run control card contain information 
required to compute the internal diameter of the flow tube from measure- 
ments taken during laminar flow runs. Each data card contains seven 
data fields each in an F10.0 format. The first 10 column field is 
Heft blank. The test section length (L) in inches is entered in the 
second data field. The elapsed time (T) in seconds required for a 
Peeactermined quantity of fluid to flow through the tube is entered in 
the third data field. The value of fluid density (SG) is entered in the 
fourth field in grams per cubic centimeters. The value of fluid vis- 
cosity (VISC) is entered in the fifth field. This value may be either 
the efflux time in seconds taken from the calibrated number 25/J680 
Cannon-Fenske viscometer, or the value of viscosity in centipoise up to 
a maximum of 5 cp. The volume (VOL) of fluid examined during the flow 
experiment expressed in milliliters is entered in the sixth data field. 
The differential test section pressure (P) which is measured in inches 
of CCL, is entered in the seventh data field. Columns 71-80 are left 
blank. 

A blank card at the end of a series will indicate that the series 


has been completed and the arithmetric average of the individually 
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computed diameters will be computed. Two blank cards at the end of a 


data deck will exit the program from the computer. Sample input data 


are located at the conclusion of the program list. 
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COMPUTER PROGRAM NUMBER 4 


Computer Program to Analyze Flow Data. 


ee 


This program was developed to analyze the experimental data taken 
during flow experiments with the variable pressure single pass flow 
system. Each series of data cards containing information taken during 
flow experiments is preceded by two run control cards. The first run 
control card contains four 10 column input fields in an F10.0 format. 

The gage number of the flow tube used during the flow experiments 
is entered in the first field (TUBE). The density in grams per cubic 
centimeters of the manometric fluid used to measure differential pressure 
is entered in the second field (SGMF). The flow diameter (D) and 
test section length (L) both in inches are entered in the third and fourth 
fields respectively. The second run control card is in a 65Al format 
on which a general description of the experiment is entered. 

There are two data cards required to record information pertinent 
to each flow experiment. The first data card contains eight 10 column 
fields in an F10.0 format. The first two fields are used to enter the 
flow tube diameter (D) and test section length (L) both in inches. The 
third field is used to enter the elapsed time (T) of the flow run in 
seconds. The value of fluid density (SG) is entered in the fourth field 
in either CGS units up to 5 grams per cubic centimeter or British 
Engineering System units for values exceeding 5 pounds per cubic foot. 
The value of fluid viscosity (VISC) is entered in the fifth field. This 
value may be either the efflux time in seconds taken from the calibrated 
number 25/J680 Cannon-Fenske viscometer, the temperature in °C of 
distilled water if used in the flow experiment between 10-99°C or the 
value of viscosity in centipoise of the fluid if less than 10 cp. The 
volume (VO) of fluid examined during the flow experiment expressed in 
milliliters is entered in the sixth data field if a volumetric measure- 
ment is taken. The differential pressure (P) is entered in the seventh 
data field. This value may be either in inches of manometric fluid when 
the VO field contains a digit, or pounds per square inch if the VO 


field is blank. The weight (W) in pounds of fluid examined during the 
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flow experiment if a weight measurement is taken, is entered in the eighth 


data field. 

The second data card contains three 10 column fields in an F10.0 
format. It is used to enter the pressure gage calibration correction 
(DEV) and maximum (PLUS) and minimum (SUB) expected error bounds 
associated with the pressure measuring instrument. 

A blank card at the end of a series of runs will indicate that 
the series has been completed. Two blank cards will exit the program 
from the computer. 

A listing of data taken during the flow experiments is located 
at the end of the program in this Appendix. Examples of output from 


this program are found in Appendices C and D. 
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CALIBRATION OF CAPILLARY TUBE GAGE 15. 


MANUFACTURERS SUPPLIED DIAMETER IS 0.054 
PLUS OR MINUS 0.2005 INCHES 

REYNOLDS COMPLTED MEASURED FLOW PRESSURE OBSERVED COMPUTED 
NUMBER VISCOSITY DENSITY RATE DIFFER. FRICTION DIAMETER 
(CP) (GM/CC) (GAL/MIN) (CIN CCL4)FACTOR (INCHES) 

80. 0.9560 0.9975 0.0 0.9 0.1992 0.0544 

eal, 0.9560 0.9975 O10 2 eS 0.9794 0.0580 

ao? . 0.9560 0.9975 0.0 So 0.0599 0.0530 
415. 0.9560 0.9975 ed 4.7 0.0385 0.052 
485, 0.9560 0.9975 abe (8 5.5 0.0330 0.0530 
Bid. 0.9560 0.9975 0.0 5.6 0.0320 0.0536 
680. 0.95460 0.9975 0.0 7.8 0.0235 0.0530 
756. 0.9540 0.9975 0.0 8.6 0.0212  “‘Oadie 
774. 0.9560 0.9975 0.0 8.6 0.9207 0.95g 
1007. 0.95460 0.9975 0.0 fess 0.0159 0.0533 

Aa 4). 0.9560 0.9975 0.0 ees 0:0148 0.0855 
ARTTHMETIC AVERAGE OF COMPUTED INTERNAL DIAMETER 0.05231 
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CALIBRATION OF CAPILLARY TUBE GAGE 13. 
HANUEAG  URERS SUPP TED PAMETIER 1S 07071 
PES SOR SNS 0 00S. Nees 


SeewOLDS COMPUTED MEASURED FLOW PRESSURE ORSERVED COMPUTED 


NUMBER VISCOSITY DENSITY RATE DIFFER. FRICTION DIAMETER 
(CP) CGM/CC) (GAL/MIN) CIN CCL4)F ACTOR (INCHES) 
136 0.9840 0,9973 Gia 0.0733 
255 0.9840 0,9973 9.0628 Diy Give 
Sere . 0.9840 0.9973 0.0442 D'. O7ame 
S42. 0.9840 0.9973 OF Oise 0.9702 
G0 . 0.9840 0,9973 0.0228 0.0718 
870 0.9840 0,9973 0.0184 0.0725 
1002. 0.9840 0,9973 0.01450 0. G7 an 
ee? 0.9840 0.9973 De One2 7 0. Oi7aes 
1396. 0.9840 0.9973 030145 0.075 
1598. 0.9840 0.9973 0.0100 0.0720 
1800. 0.9840 0.9973 0.0089 0. 7m 
1915. 0.9840 0.9973 0.9084 0.9705 
2055. 0.9840 0,9973 0.0078 0.0702 
ARITHMETIC AVERAGE OF COMPUTED INTERNAL DIAMETER 0.07158 
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CALIBRATION OF CAPILLARY TUBF GAGE 12. 
MANUFACTURERS SUPPLIED DIAMETER IS 0.085 
PLUS OR MINUS 0.005 INCHES 


REYNOLDS COMPUTED MEASURED FLOW PRESSURE OBSERVED COMPUTED 
NUMBER VISCOSITY DENSITY RATE DIFFER. FRICTION DIAMETER 
(CP) CGM7CG) GAT ae GiNece b4).F ACnoiR (INCHES) 

7 8 . 0.9828 0.9972 oe 0.9 0.0898 0.9834 
604. 0.9828 0.9972 0.0 oo 0.0265 0.9878 

a7 1. 0.98276 00,9972 0.0 Ages 0.0165 0.0865 
25 . 0.9628 0.9972 0.0 5.5 G.0131 0.0859 
1427. 0.9828 0,9972 O70 6.3 0.0112 0.0863 
1621. 0.9826 90.9972 0.0 7.3 0.0099 0.0856 
1815. 0.9828 0.9972 0.0 8.2 0.0088 0.0854 
2014. (oboe oe 10,9972 ie 9.2 0.06079 0.0852 
2149. 0.9828 0.9972 0.1 sacred 0.0074 Ono 
2154, 0.9628 00,9972 Oe sehen 0.0074 0.0844 


ARITHMETIC AVERAGE OF COMPUTED INTERNAL DIAMETER 


0.908548 
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CALIBRATION OF CAPILLARY TUBE GAGE 41. 


MANUFACTURERS SUPPLIED DIAMETER IS 0.094 
PES SOR RANG SO 0 tN RES 

REYNOLDS COMPUTED MEASURED FLOW PRESSURE OBSERVED COMPUTED 
NUMBER VISCOSITY DENSITY RATE DIFFER. FRICTION DIAMETER 
(OP) CGM/CC) (GAL/MIN) CIN CCL4)F ACTOR (INCHES) 

455, 0.9352 0.9974 0.0 1.5 0.0552 90 fees 

63S9 Oey Sec 0.9974 0.0 2+2 Ue Ooo 0.0939 

950. Oe oc 0.9974 0.0 3.4 0.0168 0.0997 

dein) 0.9352 0.9974 0.0 4.3 0.0136 0.0920 
1412 0.9352 0.9974 Cea! 51 0.01153 0.0924 
1667 0.9352 00,9974 G0 6.1 0.0096 0.9929 
1896 0.9352 0.9974 Oia 7.0 0.0084 0.0917 
2029 0.9352 0.9974 ee Te. 0.0079 0.09 
Beet iiMeliC AVERAGE OF COMe TED INTERNAL DIAMETER Oatiyoce 
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CALIBRATION OF CAPILLARY TUBE GAGE 10. 
MANUFACTURERS SUPPLIED DIAMETER IS 0.106 
PEUS OR Shins SOS 0G SSINGHES 


REYNOLDS COMPUTED MEASURED  FLOK PRESSURE OBSERVED COMPUTED 
NUMBER VISCOSITY DENSITY RATE DIFFER. FRICTION DIAMETER 
(CP) (GM/CC) (GAL/MIN) CIN COL4)FACTOR (INCHES) 

358. 0.9351 00,9973 0.0 is 0.0447 0.1084 

656. 0.9363 0.9985 0.0 1.8 0.0244 0.1063 

817. imo sc nde oos 0.0 “pel 0.0196 0.1076 
947. 0.9351 0.9973 0.0 2 0.0169 0.1064 
1164, 0.9363 0.9985 0.0 2.9 0.0138 0.1082 
1309. 0.9351 0.9973 0.0 3.6 0.0122  Oudene 
1558, 0.9351 0.9973 0.0 4.4 0.0103 0.1042 
oeB.S « GRo Soa 0.9973 Ot: 5.3 0.0087 0.1034 
2034. Gessa O57 5 ed 5.9 0.0079 0.1033 


ARITHMETIC AVERAGE OF COMPUTED 


INTERNAL DIAMETER 


07. 20520 
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DN N Sle ous F senile te 


FLOW OF TAP WATER TQ DETERMINE INTERIOP ROUGHNESS CHARACTERISTICS 
IN FLOW TURE MUMBER 45, 


FLOW TUBE DIMENSIGNS CJ]NCHES) FLOW TUBE DIAMETER 0.053% 
fest SECTION EENGT S22 


REYNOLDS COMPUTED MEASURED FLOW PRESSURE OBSERVED COMPUTED 

NUMBER VSG sil) PAVE ST 17 RATE DIFFER, FRICTION FRICTION 
CCR CEM CO) COREA) ORs 1) FACTOR FACTOR 

Oe so). 0.9762 0.9975 Oke $62.0 0.09554 0.90495 


MAX EXPECTED VALUES DUE TO UNCERTAINTY Wea 0.00558 
MIN EXPECTED VALUcS DUE TO UNCERTAINTY 256.2 0.00545 


64731. eee caee Oey eee 413.8 0.00536 0.09496 
fees xe ECTER VALUES DUE TO UNCERTAINTY 446.7 Uo 
GeeexPeCTSO VALUES DUE TO UNCERTAINTY 4065.4 C200 529 


mo. tee ye 0.9957 0,97 288.4 0.00549 0.00521 
MAX EXPECTEN VALUES DLE TO UNCERTAINTY 290.2 0200552 
MIN EXPECTED VALUES DLE TO UNCERTAINTY 285.6 0.00543 


420138. ee) aye Oe ey D7 6 Me re 0.00577 0 3viGBeZ 
fee PECTED VALUCS DUE TO UNCERTAINTY 1 oale.5 C200 ae2 
Beemer ECTED VALUES DEE TO UNCERTAINTY 167.2 0.00570 


BoA . ae 0.9927 0,54 Lye eee 0.0062nN 0.00664 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 101.4 0.00426 
MIN EXPECTED VALUES DUE TO UNCERTAINTY oe .4 0.09606 

eo? . Dee 2 0.9957 Ona 2 B52 0.00663 0.00 ec 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 67.5 0.0649 
MIN EXPECTED VALUES DLE TO UNCERTAINTY 65.6 O70 St 

foo 04. Pee au 60,9967 0.36 pie C.00725 OQ,UGge.e¢ 
MAX EXPECTER vALUES DLE TO UNCERTALYTY 51.9 O ti 22 
MIN EXPECTED VALUFS DLE TO UNCESTAINTY see 0.00709 





FLOW OF PURE JF4 TO DETERMINE INTERIOR ROUGHNE 
IN FLOW TUBE NUMGER 45, 
FLOW TUBE DIMENSIONS (INCHES) FLOW TURE DIAMETER 00,0533 
TEST SECTION LENGTH 
REYNOLDS COMPUTED MEASURED FLOW PRESSURE 9RSERVED 
NUMBER VISCOSITY DENSITY RATE NIFFE, FRICTION 
(C°) Gi Cie? ay) CPS) Pn Gaim 

64257 fe 7 46 0.7644 0.73 128.14 G.-C 5> 
Meee XPECTER VALUES DUE TG UNCERTAINTY 129.1 0.90569 
MIN EXPECTER VALUES DUE TQ UNCERTAINTY ee ico 0.090554 

B77 55, 0.7438 00,7644 0,62 OF = 7 0.00599 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 98.4 OOD ee 
PEMPGXPECTER VALUES CUE TC UNVERTSINTY 95.8 Dee Sas. 

99444, Geass Ss 0.7654 0,48 63.9 0.09626 
MEX EXPECTED VALUES DUE TC UNCERTAINTY ou GC .99537 
Pewee XPECTED VALUES NUE TO UNCERTAINTY Byieo 0.) Sie 

29270. 0.7942 OLB oO? rear! 69.1 G0 0s14 
MAX EXPECTED VaLUES QUE TO UNCERTAITY yon 0.99419 
MIN EXPECTED VALUES DUE TO UNCERTAINTY As «2 0 eS 32 

25489, Oe7 942 eo, Omoe 55, 0.99449 
Pee x PECTED VALUES DUE TQ UNCERTAIATY 54.4 0.00656 
BSeermee Xe eCTED YALVYUES DUE TO UNCERTAINTY 64,0 0.90430 

P2iis6. 0.7942 0.7657 0,39 41.4 0.90441 
Mex EXESECTEN VALileS DUE TO UNTERTAINTY ZAG 9.10548 
MIN EXPECTED VALJES DUE TO UNCERTAINTY 40,2 i, Coase? 

foac). 0.7942 Ron? oe 7 OF ees Olas 0.00635 
Peeee xP =CTER VALUES DUS TO UNCERTAINTY 220 0.00542 
MIN EXPECTED vALUES DUS TO UNCERTAINTY 2152 Oo Grey, 


ANALYSIS OF FLOW 


DATA 


0. 


0 


0 


Be 


0 


0 


0 


oan 


SS CHAP ACTERTS Ties 


5.405 


COM eo 
FoR TiC ean 
FACTOR 


ee we Oe | | 


090545 


»09604 


J169%3 


~99626 


~09449 


007 G2 
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ANAS lon Ur io Oul a A tn 


FLOW OF PURE J®4 TO DETERMING INTERIOR ROUGHNESS CHARACTERISTICS 
[IN FLOA TURE NUMBER 415, 


moO TURE DIMENSTONS © | NiGieS) FLOW JUBE DIAMETER 20,0533 
TEST SECTION LENGTH 5,405 


REYNOLDS COMPUTED MEASURED FLOW PRESSIRE ASSERVED COMPUTED 
NUMBER ViisG@misi ty DENSt TY RATE HlFPEQ. FRICTION. FRIIS. 
. CCa ir (GCM /OC}. (Sky ee tes) FACTOR FACTOR 

Rep O7 . G27295 eae 0.64 Ae 0.005580 0 70seae 


Beeeme PECTED VALUES DUE TCO UNCERTAINTY su) Bic 0.00568 
Meme =eCTED VALUES DUE JO UNCERTAINTY Aes 0. C0 >a 


oe . V7 eae Oro rear J 7 Ane 9.99547 0.99590 
Seema se -CTED VALUES DUE TO UNCERTAINTY 17227? 0.0)554 
Pee mee OTED VALUES DUE TO UNCERTAINTY 7 a 2 G2 


i Soysiem O27 ae Oy 2640S ee 314.6 0.90519 0.09480 
Pemex PECTED VALUES DUE TO UNCERTAINTY Sie 0.99524 
fee PECTED VALUES DUE TO UNCERTAINTY Si 0.90512 


SSO SS . De Oi os ei a5 72 0.00520 0.00 
femme PECTED VALUES DUE TO UNCERTAINTY 439.9 0.00523 
Beemer eCTED VALUES DUE TO UNCERTAINTY 432,4 0.00514 


e2ce1. Og OS O27 6<4 Rta) 2 247,4 6.99547 Oeaae ie 
Siem EOTED VALUES DUE TO UNCERTAINTY 250.4 0,99554 
MIN EXPECTED VALUES DUE TO UNCERTAINTY Pe Aad 0,9954n 


poo? Oo 0.7644 0,88 1eSed., 5 0,99552 O2e3570 
Pemex PECTED VALUES DUE TQ UNCERTAINTY 1:8 Se. 70S 64 
MIN EXPECTED VALUES OUS TO UNCERTAINTY Pate 9,90549 
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ANALYSJS OF FLOW DATA 


FLOW OF TAP WATER TO DETERMINE INTERIOR ROUGHNESS CHARACTERISTICS 
INF POPS UBE  MUMEE Ras: 


FLOW TURE DIMENSIONS C]NCRES) FLOW TUBE DIAMETER 0.0716 
TEST SSECTION EENGTH 7a ee 


SeeNOLDS COMPUTED MEASURED PLOW PRESSURE OSSERVED COMPUTED 

NUMBER VISCOSITY DENSITY RATE DIFFER, FRICTION FRICTION 
Coe) GCu7 ce. (GAL/SMIN) (PST) FACTOR FACTOR 

9665. ete 64 00,9987 OG 424,8& G0,090469 0.00449 


MAX EXPECTED VALUES ODUE TO UNCERTAINTY 437.4 ee 2 
MIN EXPECTED VALUES OUE TG UNCERTAINTY 43n,9 0,00465 


83247, 1.0691 00,9987 Zo 2 Se Oia 0.00479 UeOtece 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 3326 C,00482 
MIN EXPECTED VALUES DUE TO UNCERTAINTY Sa 0 i 00474 


yeo68 . i) Oo 1 De ie? Wea’ Coe a8 0.00488 0.00481 
Meee xXPECTED VALUES DUE TO UNCERTAINTY A 0.00494 
Pee XP=CTED VALUES DUE TO UNCERTAIMTY 25a 7 0.00485 


eo? . 0.95564 0.99689 1,39 6c 7, 0.00500 OQ 29GSGR 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 162.2 9,00505 
MIN EXPECTED VALUES DUE TO UNCERTAINTY 169.6 C.00494 





ANALYSIS OF FILOW 


DATA 


94 


Te Oe 


COMPIJTED 
FRICTION 
FACTOR 


— ee oe oe 


0.00445 


3.09478 


0.09528 


). yom 


0.0359 


0.09595 


0. 000 


) Jae 


FLOW OF PURE JP4 TO DETERMINE INTERINR ROUGHNESS CHARACTERISTICS 
IN FLOW TUBE NUIRER 42, 
Pio TURE DIMENSIONS (INCHES) FLOW TURE DIAMETER 0.07164 
TEST SECTION LENGTH 
REYNOLDS COMPUTED MEASURED FLOW PRESSURE ABSERVED 
NUMBER VOESeaicnliny suens 1 1y RAT= DIFFER, FRICTION 
(CP) CGEr7CO) (CG le eee ST) FACTOR 
99800. 0.7124 O.. 7 SAD Zl 269.3 0.094665 
Pee EXPECTED VALUES DUS FO UNCERTAINTY 2756 On00Mizg 
MIN EXPECTEN VALUES DUE TO UNCERTAINTY 264.7 0.00441 
74795. 0.7124 0.7410 4.59 159.4 0.900499 
MAX EXPECTED VALUES OUE TO UNCERTAINTY 140.7 0,99495 
MIN EXPECTED VALUES DUE TO UNCERTAINTY 15520 0,90489 
50394. O07 baa 0,74619 eae ee 0.0N522% 
MAX EXPECTED VALUES DUE TO UNCERTAINTY ae 0.00529 
MIN EXPECTED VALUES DUE TO UNCERTAINTY 75.9 0.005909 
55823, Oa 7e4s 9.7644 jee ao 0,0051¢ 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 198.4 0200523 
MIN EXPECTED VALUES DVE TO UNCERTAINTY te eet Oe OM Sar 
462670, 0.74546 0,74544 9,97 67.0 0,0955% 
MAX EXPECTED VALUFS DUE TO UNCERTAINTY eo 0,90540 
MIN EXPECTED VALUES DUE TS UNCH2TAINTY 65.0 0.00536 
Buea? . 0.74546 0.7444 ee S76 hw: 
MAX EYXY2ECTEN VALUES DUE TO UNTERTAINTY So 0.00589 
MIN EYPFCTED VALUES QUE TO UNCERTAINTY 345.1 0.00557 
S7016. 0.7544 06,7544 0,384 51.5 0.09545 
Pee xPSCTeDR VALUES OLE TO UNCERTAINTY Bo 9.9957% 
MIN EXPECTED VALUES DUE TO JINZERTAINTY 59.9 0.00548 
Be251. 0.746464 0,7444 0,59 a 0.00544 
MEY EXPECT=EL VALUTS DUS TI UNCERTAINTY 94.2 0.99573 
MIN EXPECTER VALUES DUE TO JNCERTAINTY 2 ae 0200 550 
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ENN So OE Sea aay A 
FLOW GF TAP WATER TC DETERMINE INTERIGQR ROUGHNESS CHARACTERISTICS 
} IN FLOW TURE NUMBER 10. 
FLOW TURE DIMENSIONS CYJNCHES) Pov e= The DTAME imeR Oe 
TES! SECTION LENGTH a0y oe 
NUMBER VVSG@sct TY DENSITY RATE DIFFES, FRICTION FRI Giga 
(C2) Cec) (GALSMIN) (PS]) FACTOR FACTOR 
yoy 1S . ieee 2 90,9975 A Onl) AGS 0.00578 £40.04 77 
MAX EXPECTED VALUES NUE TO UNCERTAIPTY 2 SS 0.00522 
MIN EXPECTED VALUES DUE TO UNCERTAINTY 120737 O00 507 
56710. dee 0 leer 00,9975 aos Von) 0.00570 Jed S es 
MAX EXPECTED VALUES NUE TO UNCERTAINTY Loe 60574 
Pepe XPECTED VALUES DUE TO UNCERTAINTY pared ) -00S5 7 
a5G//. aloe) ee 0.9975 ion S252 0.0061n OU2voSa 
MAX EXPECTED VALUES DUE TG UNCERTAINTY 52.9 0, Goes 
feeeemerxPCCTED VALUES DUE TO UNCERTAINTY S41.1 0.00593 
noe 55, Ie Oe Spe 00,9975 4,12 Staeno 0.990515 U0 aie 
MAX EXPECTED VALUES DUE TO UNCERTAINTY ohare G0 a5 
MIN EXPECTED VALUES DUS TQ UNCERTAINTY So, 7 0.00 S129 
84734, 1.05384 0,997 Cerne 17 te > 60,99527 0.99464 
MAX EXPECTED VALUES Dus an “UNCERTAINTY is/ dag OG os oo 
MIN EXPECTED VALUES OUS TC UNCERTAIS: TY 148.2 6.00320 





ANAISY S loa 0b s Pigoil 


DATA 


96 


COMPUTED 
FRI Gl G 
FACTOR 


0.00456 


0.00480 


OO Saas 


0.00540 


FLOW OF PURE JP4 TO DETERMINE INTERIOR ROUGRENESS CHARACTERISTICS 
IN FLOW TUBE NUMBER 49, 
TEST SECTION LENGTH 10.609 
ReyNOLDS COMPUTED MEASURED FLOW PRESSURE OBSERVED 
NUMBER VISCOSITY DENSITY RATE DIFFER, PRC Gs 
(CO) CG Ge) (GAL AMIS (PS]) FACTOR 

Poo. Gme7i47.5 0.7633 2.98 128.9 eo OSs 
Meee APECTCOD VALUES DUE TQ UNCERTAINTY 30004 OCS s2 
MIN EXPECTED VALUFS DUE TC UNCERTAINTY 1 iG ee CoS 2h 

oO. 0.7475 (276355 2.42 Oi ao 0; 20S43 
Mex EXPECTED VALUES DUE TO UNCERTAINTY §6.3 0,860555 
MIN EXPECTED VALUES PLS TO UNCERTAINTY &5.0 CPSs 4 

BIoo oi. Oe 7647 G,74Sa 8 63.9 0.006450 
MAX EXPECTED vALUES DUE TO UNCERTAINTY 64,5 0,00459 
MIN EXPECTED VALUES DUE TO UNCERTAINTY 62.0 CC. 00430 

aot 1.5. Gn 7ode, 0,76S8 5G 38.3 0.090589 
MAX EXPECTED VALUES DUE TO UNCEPTAINTY Zk, / O40 oS 
MIN EXPECTED VALUES CUE TO UNCERTAINTY S74 i C2005 75 

64065. 0.7647 0.7638 Totee 20.23 0.0052? 
MAX EXPECTED VALUES DUE TO UNCERTAINTY ote C.00588 
MIN EXPECTER VALUES DUE TO UNCERTAINTY Pile 0 U0 S45 
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Ae Yoel eS ee ames @imaa 1k 


LAMINAR FLOW TEST USING PURE TAP WATER 
IN FLOW TUBE NUMBER 43, 


FLOW TUBE DIMENSIONS CINCHES) FLOW TURE DI] AMETER 0.0716 
TEST SECTION LENGTH 7.109 


REYNOLDS COMPUTED MEASURED FLOW PRESSURE ORSERVED COMSUTED 
NUMBER ViES@aSITY DENSITY RATE DIFFER. FRICTION FRICTION 
(CP) (GM/CC) (GAL/MIN) ¢PST) FACTOR FACTOR 

aD, 6.9840 0.9973 One Oa20 0.04298 0.94418 

Bt?! . 0.9840 0.997% Oot eae 0.03443 0.03185 
mS), 0.9840 be oor Wye Brea! 0,02245 Oi d22e. 

884 0.9840 Cc fas 202 ie 0.01730 9 foaretee 
1002. 0.9840 00,9973 G2 et 0.01598 0.9907 
eee 0.9849 0.9973 OOS, Ger2 0.01298 9.01278 
imo? >. 0.9840 USO es OS 0.2 0.01151 Uy Gi 
mou . 0.9840 6.9973 0.04 0.2 0.00974 CC. o7CCe 

fie 6. 0.9846 O59 97S 0,04 0.5 0.00947 @.90900 
1884, 0.92849 00,9973 0.04 Os 0.00°05 9.90849 


5/ 0.99722 


wy 


ao i6 . 0.9840 De 9 oes ea N.S 0.90 





FLOM 


REYNOLDS 


NUMBER 


LAMINAR FLOW 
LN Sie 


TURE DIMENSIONS 


COMPUTED 
VISCOSITY DERSITY 
Cel 7 CCD 


ae FF oe OP ee —_ “ae wow | ae 


0.9566 
eo oet 
Um OU 
0.9560 
(eo 208 
Or go SU 
0.9560 
Og 556 


G2 60 


ANE Y SiS 


MEASURED 


Ones 
O. 9975 
0.9975 
0.9975 


G.9975 


(]NCHES) 


OF FLOW 


TEST USING 
eee 


FLOW 


Fie hiw 
QATE 
(GAL/SMIN) 


Ce — i eel 


DATA 


Bes SU 


OOP ee 
(eo) 


PA ke Ae 
NUM Ee Ro 41 5), 


LENGTH 


Cn Seu =e 
PR iis 
pee rati 


Or oye 
Ce oa 
0.024285 
CSS ie 
G2 0St2e 
0.020958 
Gt Ao 
Oe Ges es 


GeO ee 


98 


TUBS DIAM ete 0 oS. 


ieee © 1 LON 5,406 


COMPUTED 
FRICTION 
FACTOR 


0.96029 
0.03864 
0.02568 
US Sas 
0.05182 
0.32063 
0. U2n27 
0.0 i4ea 


03a 





ag 
APPENDIX D 


GOreC ek ey kt hes OU 
RESULT Seer FE LOASDATA. TAKEN 
DURING DRA Ge eDUuC TT Toner xAPERIMENTS 


(FLOW OF 50 wPPM PIB IN MILITARY JET FUEL) 





FLOW 


REYNOLDS 
NUMBER 


80079. 


100 


ASAE SPS. 2 0iF =F Gea hay ok 


FUON OF JE) Ue te ee eC oc e eeeeeCR Oia PASS) 
IN FLOW TUGE NUMBER 414, 


TUCE SU UMENS TONS [Gee Si ee eee | ie ER 0s, (ins 
Test isie@ 7 LON EENGT4 5,70 


COMPUTED MEASURED FLOW PRESSURE NASERVED COMPUTED 
ViSm@eliyY DENSITY RATS DIFFE®, FRICTION FRiGiiee 
CGR j= (EM/CTC) (GRie/ Ish) (PST) FACTOR FACTGR 

0.7739 C7 Ges 1.36 on 6 9.00399 0.00479 


MAX EXPECTED vALUES DUE TO UNCERTAINTY 323.4 0.00¢04 
femeesPECTED VALUES DUE TO UNCERTAINTY 3144.4 0.00395 


56276. 


0.7739 0,7685 0,96 LS Sh 0.00295 0.09514 


Deemer? ECTED VALUES OLE TO UNCERTAINTY Des G Corker 
Reeve x? CTE) VALUES DUE TO UNCERTAINTY WS ok OS GO Sae/ 


43909. 


Die JS 00,7685 07S 9°,.2 0,0041> 0 .eieee 


Meee xPeCTEN VALUES DUE TO UNCERTAINTY 199.2 C,904264 
mieex-=CTED VALUES DUE TO UNCERTAINTY 97.9 0.0040¢ 





OE 


ANALYSTS OF FLOW DATA 


FLOW OF JET FURL WITH DEGRADED PIR (CTHIFD PASS) 
IN FLOW TUS= NUMRER 46, 


a0. TURES DINENSIUNS () VOne EG See ieee Or 
TEST SS=SCILOM BENGi OS Sie 


BermNoOlaS COMPUTED MEASURED FLOW PRESSURE @ESERVED COMPRiiiay 
NUMBER Veeslem@ienay y “ENS ] 7 Y Ae DTFREr. FRICTION Riv maa 
Go OGM 7 GCC) GG Aly te yee FACTOR FACTOR 


78847, 0.7654 0,74672 ts QR1,4 0.00370 0.09472 
MAX EXPECTED VALUES DUE TO UNCERTAINTY 2842.8 Ge Oe a4 
MIN EXPECTED VALUES DLE TO UNCERTAINTY DE eh 0.007464 


erg S . 0.7454 C17 ae 2 ese 0.00256 0.00595 
MAX EXPECTED VALUES DUE TO UNCERTAINTY aie 0.00379 
MIN EXPECTED VALUES DUE TO UNCERTAINTY Wet. 0 0.00359 


47755. 0.7654 Ca Se eco ee 0.00363 0. 00545 
MAX EXPECTED VALUES DUE TO UNCERTAINTY ac 6.00347 
MIN EXPECTED VALUES DUE TG UNCERTAINTY 9¢.90 0,00855 





FLOW 


REYNOLDS 
NUMBER 


8479. 
12084. 
14941. 
ieee 5 . 


20212. 


N 
Nm 
n 
© 
Nm 
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AN REO OF EIeG hb ien 


FC WecQ © su'6 lei 1) east 
J 


PIR (SECOND PASS) 
Nae Woes - 


TURE DIMENSIONS CINCORES) FLOW TIRE DIAMETER 0.050 
Test SeCriON, LENGTH 5. cio 


COMPUTED MEASURED FLOWN PRESSURE OBSERVED COMPUTED 
VISCOSITY DENSITY RATE DIFFEP., FRICTION FRICTION 
(G2) 5 (GM/CC) (GAE7MIN) ~(PS1) SIG hes FACTOR 
0.75581 0.7869 0.14 5.6 0.00647 0.90824 
fe7o3i 0.7569 te o.i - 0,00527 (8 iia 
Gg >3e) 0, 7669 es 4267 0.00473 “ONGizae 
0.7581 0.7669 0.30 46.1 0.00419 0.90684 
0.7561 43,7669 0,34 4€.5 0.00377 0.99643 


Or ae 0.7669 Oey 22.7 0.00377 0.90647 
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ANALYSIS OF FLOW DATA 


FLOW OF SOQWP2M PIQ JN JET FUELCFIRST PASS THRU ORIGINAL SYSTEM} 
IN FLOW TUBE NUMBER 45, 


meow [URE DIMENSIONS CINCHEG)  SAEOn TUS DPAMWetes 0 Uiess 
he Sieoce) | CN Mee NG eo ae 


meme oes COMPUTED MEASURED Pletole. PPESSURE GBSERVED COMPUTED 
NUMBER Vee Te ets | TY RATE DIFFER. FRICTION  —FRiUWean 
(CP). (GM/CC) (GAL/MIN) (CPST) FACTOR FACTOR 
Zoo 0. Oe o- oe 0,48 19.4 0.09240 0.90624 
om 1 4 . Oy cnea 0,7652 cre a are, 0.00784 Uv. odoa 
ooo 1. te? oe2 OO Ora B.6 0.00360 0.00766 
a) >. OV Oe2 0 7652 Ore as 0,00SS6. 0.90723 
go! 9 . Oi) elie Bf 52a2 0,48 Dic. 9 0,00224 0.90609 





FLOW OF SOWPOM PI8 


BEOA TURE DIME ySICNs 


BeyNOLOS COMPUTED 
NUMBER 
CEP) 


55005. M7705 
MAX EXDECTED VALUES 
MIN EX®5CTED VALUES 


42242, lie 0 
MAX EXPECTED VALUES 
MIN EXPECTED VALUES 


35097. Ora 7 G2 
PeeeeexPECTED VALUES 
femme xe = CTcD VALUES 


oo 17 3. Diane 
femmes? ECTS) VALUES 
pyevex@=CTSD VALUES 


eoea3 . Ura > 
Pepe eCTS) VALUES 
Pee XP ECT&OND VALIES 


ViesGoeai VelNSi ty 


ANALYSIS OF FLOW DATA 


NS Se Ue ee Tee cs 
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THRU ORIGINAL SYSTEM) 
IN FLOW TUBE NUMRER 45, 


9.406 


COMPUTED 
FRICTION 
FACTOR 


O70Cs 52 


0.00587 


0.00600 


dv tees 


(INCHES) FLOW TUBE DIAMETER 0.0533 

TEST SECTION LENGTH 
MEASURED FLOM! PRESSURE NABSERVED 
RATE DIFFES, FRICTION 

COMA) (GAL/MIN) (PS]) PA hia 
0.27670 0,90 eG OnieesS 
QUE TO UNCERTAINTY eG 1.60227 
NLE TG UNCERTAINTY 74.0 0.09218 
09,7679 OP Af S 0.09182 
QUE TO UNCERTAINTY ac. 9 0.00185 
QUE TO UNCERTAINTY 3°,0 0.99176 
07570 0,56 95.9 Gave als3) 
QUE TO UNCERTAINTY 26.2 0.09494 
DUE TO UNGER Tet Ty 25.0 0.09184 
6.7470 Casa oe 0.09184 
NU= TO UNCEFPTAINTY Oat 0.99184 
NUE TO UNCERTAINTY Dee 0.00177 
0.76/79 6.40 heal 0.00189 
DUE TO UNCERTEINTY tices Or Oo 
DUE TO UNCERTAINTY 12.9 G.0018% 





KivALY ols. OF leo 


DATA 


105 


7.109 


COMPUTED 
FRICTION 
FACTOR 


— = ae TS ae oe oe 


0.00411 


OG Oia 


0.990461 


0.0047 


OO Sag 


FLOW OF SOWPOM PIB IN JET FUEL CFIRST SASS THQIU ORIGINAL SYSTEM) 
IN FLOW TUBE NUMRER 13, 
FLOW TURE DIMENSIONS CINRChKES) FLOW TURS DIAMETER 9.0716 
TEST SECTION LENGTH 
ReereOLD)S COMPUTcCS ic ASURED FLOW PRESSURE ORSERVED 
NUMBER VISCRST TY <UeNSITY RATE DIFFER, FRICTIIN 
Cola Col 7 CC ) SaRGGsie ia SG eames cin) Factor 

13/7699, O27 761 Oe he 5 Sl 6 S742 0,0900268 
MAX EXPECTED VALUES DUE TO UNCERTAINTY $74.9 On 2 90 
MIN EXPECTED VALUFS QUE TC UNCERTAINTY S70 0.00285 
iol 7 . Hie 6 li ie Ses 2 Od Donia? Omi? > 4 
MAX EXPECTED VALUES DUE TO UNCERTAIATY Bo 44 ae 5 4 
MIN EXP®ECTER VALUES DUE TO UNCERTAIATY Dot. Uae oa 
84929. 047237 0.7655 Be) ek 47 noice 
MAX EXPECTELR VALUES DUE TO UNCERTAINTY 19s, c 0,09224 
MIN EX®SCTEN VALUES DUE TO UNCERTAINTY eee Opera 7 
e744). Gees 7 0.7655 joy 62.9 G,001 9s 
MAX EX°@ECTEN VALUES DUR TQ UNCERTAINTY 63.7 0.99198 
MIN EXP=CTED VALUES QUE TG UNCERTAIATY 614) 0.60199 
54703. 0.7837 0.7455 line ore Tid foe 
Meee X@ECTEDN VALUES NUE TO UNCEPTAINTY 5.0 Odo 
MIN EXPECTED VALUES OUE TO UNCERTAINTY Sey 9.09175 





106 


ANALYSTS OF FLOW DATA 


FLOW OF SOWPPM P13 IN JET FUEL CFIRST 2®4SS THRU ORIGINAL SYSTEM) 
IN FLOW TUBE NUMBER 15. 


meow. TURE DIMENSIONS CIMCRES) “FLOW TUB= DIAMETER 0.0558 
Teo) See, lON LENGE “>i4 ue 


REYNOLDS COMPUTED MEASURED FLOM PRESSURE S8SERVED COMPUTED 

NUMBER ViES@mciryY LOENSITY RATE DiBPe®. SRICTION “FRiGian 
(C2) |. (6BM/CC) (GAL/MIN) ¢PS]) FACTOR FACTOR 

eo 7 2. 0.72366 0.7633 ea Ay en 0.90348 0.00435 


MAX EXSECTED VALUES NUE TO UNCERTAIATY 473.9 0.00350 
MIN EXPECTEN VALUES DUE TO UNCERTAINTY 465.9 0.00344 


Oe 39 . 0.7368& 17665 41-1, 50 286.4 0.090295 0.00454 
MAX EXSECTED VALUES DUE TO UNCERTAINTY 287. 0.09299 
MIN EXPECTED VALUES DUE TQ UNCERTAINTY 283%,6 P2009 297 


pease. ore G7 65S fee 17 04 0.00287 0.00481 
MAX EXPFCTED VALUES DUE TO UNCERTAINTY a eee 0.04290 
MIN EXPECTED VALUES DUE TC UNCERTAINTY ie ee 0.99281 


56841. Giese soc 9.7638 0,93 To ae 0.99201 0. U0 oir 
MAX EXPECTEN VALUES DUE TC UNCERTAINTY 74.9 oO, O02 04 
MIN EX®=CTED VALUES DUE TO UNCERTAINTY 7 oat) 0.091964 
S00 0.7566 1760s Cano fied 0.06749 0205635 


1352. Oey S57 uf oe S G2 0.3 G.01SSS 5 0 oe 
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ANALYSIS OF FLOK DATA 


FLOW OF JPe WITH DEGRADED PIR (SECOND PASS THRU MODIFIED SYSTEM) 
IN FLOW TUBE NUMBER 45. 


mma! TUBE DIMENSIONS C]NCMES)) GOW TUR] DIAMETER 0.0233 
(ese ON CRENG Ts. 2 Uio 


Meee OlDS COMPUTED MEASURED PON DEESSURE OBSERVED GOM=UieD 

NUMBER VISCOSITY DENSITY RATE Niehes. “FRICTION, “a Cian 
(C>) _ (GM/CC) (GAL/MIN) ¢PS7) FACTOR FACTOR 

eyo. . Ou Oe GO 7 G9 Oe 2 23.2 0.09655 2790GCu 

14744, OR ok 0,7669 ee 13.2 0.00698 Us 07a 

120350. oo 0,7669 0.20 13.4 0.080751 0 0nviS 

eoa7, Oo 6 a 0,7669 O55 7.9 0.00826 > 0 l[Uiieae 


eeewO . O72 3) 00,7669 Cee 13.4 0,00774 0.0%736 





Piso OF 


FLOW 


REYNOLDS 
NUMBER 


Ce ee ee 


S50WPPM P18 


LUBE DivieNs TONS 


COMPUTED MEASURED 

VS BENS Ty 
CCE) (GM/CC) 
0.7624 1.7654 
Jie? feo li, 766 
Ooo 0,7654 
0.7624 0,7654 
C7 0.7654 


TN 


ANALY Se Oh ak 0 


IN FLOW 


este 


(TNCFES) 


PSO ee 


DATA 


Teste seo) ON WeeNGir 


FLOV' 
RATE 
(GAL/MIN) 


oe Ge em 


Cod 


PRESSURE OSSERVED 
DIFFER, 
Gea) 


FRICTION 
FACTOR 


-_— Swe ew eee ee 


6.00477 
0.00487 
0.00474 
0.00542 


GO. U0sat 
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JET FUELCFIRST PASS THRU MODIFIED S¥SGua yD 
NUM Bena. 


PAM EIER 020245 
54 405 


COMPU s 
PREC Rion 
FACTOAR 


eee Te SY ee 


0.00678 
0.0 0m7 2 
0 wise > 
0.00702 


0 0M as 





mii OF 


FLOW 


REYNOLDS 
NUMBER 


2270. 
1491. 


a7 5 . 
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KN AY Sues OF oF SO Peieiae 8 


GONPEM PI3 IN JET FUEL CFIRST PASS THRU MODIFIED SYSTEM) 
IN FLOW TUBE NUMRER 45, 


TURE DIMENSIONS CITRCKES) FLOW TURE NIAMETER 9.0538 


TEST SSeG TON LENG Te 5.405 
COMPUTED MEASURED FLOW PRESSURE ORSERVED COMPUTED 
VaECCOStal vate Ss | 7 Y RATE DIFFER, FRICTION FRICTION 
(CP) (GM/CC) (GAL/MIN) (€PS]7) FACTOR FACTOR 
a7 578 0.7642 0.23 Le ac 0.90498 0.09692 
yoy 2 0.7642 0,19 duet 0.90634 0.99764 
Uy Coe Oeed > 6.8 G@,09736 Ciooede 
0.7578 ey re pe 4.8 0,00824 0.09866 
D577 8 00,7642 ao Ze 0,00851 0.33950 
tiny 276 oa O77 1.9 O,00943 “v.doeer 
07576 0,74642 OF G6 1.3 0.00990 UU. 0222 
0.7578 0.7642 one N.9 1.00940 0.0 tdee 
Uae G 0.7642 C.04 c.8 0.00889 0.011902 
Q@. 7576 00,7642 C.04 N.S 0.00841 0.011443 
Ooi s 9.7442 ene Dans O02 546 0 oa 
7 37 0.7642 Sean ieee 0.01954 0.91829 














oj 7 % 
PE ey Pa 
ce a Oe 
RT 7} Wi 
' g rei 


an 


3 
Hy 


oa if 


Studies to increase the flow rate of mil 


| HNN HIT Hit i Wi INIT Will 
AHA AER En ill 
i i! WAL ML AUT} WLU lt | Mt 


3 2768 001 94773 2 
DUDLEY KNOX LIBRARY 


ree 


Py; 








